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ABSTRACT  OF  THESIS 

CALIBRATION  OF  THE  SCANNING  SPECTRAL  POLARIMETER  AND 
MEASUREMENT  OF  THE  SKY  LIGHT  POLARIZATION 


The  Scanning  Spectral  Polarimeter  (SSP)  measures  6  optical  properties:  flux,  unpolarized 
radiance,  and  four  polarized  radiances  (parallel,  perpendicular,  right  hand  circular,  and  left 
hand  circular)  through  the  spectral  region  from  400  nm  -  4000nm.  The  SSP  was  designed 
for  the  measurement  of  reflected  sunlight  and  a  complete  description  is  provided.  This 
thesis  asses  the  ability  of  the  SSP  to  measure  the  clear  sky  spectral  polarization.  Initial 
calibration  procedures  and  results  are  discussed.  Instrument  characteristics  are  provided  to 
include,  field  of  view,  transmission  function  for  the  diffuse  channel,  and  detector  calibration 
coefficients. 

A  plane  polarized  radiative  transfer  model  is  used  to  study  the  effects  local  conditions 
have  on  the  sky  light  polarization.  The  impacts  of  changing  Rayleigh  optical  depth,  surface 
albedo,  solar  position,  and  haze  loading  are  examined.  The  SSP  is  used  to  measure  the 
spectral  sky  light  polarization  and  the  results  are  compared  to  those  values  predicted  by 
the  model.  The  SSP  can  determine  sky  light  polarization  to  within  10%  error  with  respect 
to  the  model  predictions  and  was  capable  of  resolving  the  effects  aerosol  scattering  and 
surface  reflection  have  on  sky  polarization. 

Vincent  T.  Ries 

Department  of  Atmospheric  Science 
Colorado  State  University 
Fort  Collins,  Colorado  80523 
Fall  1995 
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Chapter  1 


INTRODUCTION 


1.1  Motivation 

The  earth’s  climate  is  controlled  by  the  spatial  and  temporal  variations  of  radiation 
received  from  the  sun.  Incoming  solar  radiation  is  balanced  by  energy  that  is  reflected 
and  emitted  from  the  earth-atmosphere  system.  Variations  in  outgoing  radiation  are 
predominately  influenced  by  clouds,  whose  temporal  and  spatial  variations  are  not  well 
understood.  Even  less  understood  are  the  radiative  transfer  properties  of  clouds  (Takano 
and  Liou,  1989)  and  how  these  properties  relate  to  cloud  microphysics  and  ultimately  to 
radiative  heating  and  cooling  (Tsay  et  al.,  1994). 

Because  of  the  complexities  in  modeling  the  radiation/cloud  interaction  there  is  a 
requirement  for  high  quality,  high  resolution,  spectral  observations  of  cloud  and  aerosol 
properties  which  can  be  used  to  determine  areas  where  our  understanding  is  lacking.  It 
is  also  important  to  compare  albedos  at  both  visible  and  selected  near  infrared  windows 
Eis  this  relationship  provides  microphysical  information  including  particle  size  (Nakajima 
and  King,  1990)  and  (Stackhouse  et  al.,  1994). 

Sky  light  intensity  and  polarization  calculations  were  pioneered  by  Chandrasekhar 
and  Ebert  (1954).  The  first  studies  involved  only  Rayleigh  scattering.  As  computer  power 
increased,  methods  were  developed  that  included  scattering  from  large  spherical  particles 
(Dave,  1970),  then  haze  and  aerosols  (Kattawar  et  al.,  1976),  (Hitzfelder  et  al.,  1976), 
and  (Dave,  1978). 

Coulson  (1977)  used  skylight  measurements  to  determine  atmospheric  turbidity.  The 
effects  on  sky  intensity  and  polarization  from  haze,  ice  crystal  precipitation,  and  volcanic 
clouds  have  been  studied  by  Deliver  (1988),  Fitch  and  Coulson  (1983),  and  Coulson 
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(1983).  The  polarization  and  intensity  of  sky  light  from  the  zenith  sky  has  been  used  to 
determine  upper  tropospheric  and  stratospheric  turbid  layers  (Coulson,  1980),  (Coulson, 
1981),  and  (Belying  and  Lu,  1988). 

Interest  in  sky  light  intensity  and  polarization  has  again  peaked  with  the  advent  of 
the  Polarization  and  Directionality  of  Earth’s  Reflectances  (POLDER)  experiment  that 
will  fly  on  the  Advanced  Earth  Observing  System  (ADEOS)  satellite,  scheduled  to  launch 
in  1996  (Deschamps  et  al.,  1994).  In  preparation  for  this  experiment  an  airborne  version 
of  the  POLDER  began  flying  in  1990.  Data  from  the  airborne  version  of  POLDER  have 
been  used  to  estimate  aerosol  loading  (Deuze  et  al.,  1993),  determine  liquid  phase  of  cloud 
particles,  and  derive  cloud  top  altitude  (Goloub  et  al.,  1994). 

1.2  Why  the  Scanning  Spectral  Polarimeter? 

An  instrument  was  developed  that  is  capable  of  measuring  the  spectral  properties  of 
solar  and  near-IR  radiation  reflected  from  clouds  at  a  resolution  fine  enough  to  adequately 
describe  the  spatial  variation  of  this  energy  reflectance.  The  Scanning  Spectral  Polarimeter 
(SSP)  provides  radiance  measurements  as  a  function  of  wavelength  from  400  nm  through 
4000  nm.  As  presently  configured  the  SSP  provides  spectral  measurements  which  are 
unique  to  the  remote  sensing  of  cloud,  aerosol,  and  land  surfaces  in  the  solar  wavelengths. 
From  these  data,  the  SSP  will  provide  a  wealth  of  information.  The  spectral  nature  of  the 
data  provides  information  on  cloud  optical  depth  and  effective  sizes  of  cloud  water  and 
ice  particles  (Stackhouse  and  Stephens,  1991).  The  polarization  measurements  are  useful 
in  determining  cloud  optical  depth  and  microphysical  properties.  From  the  reflected  and 
transmitted  radiances,  the  fluxes  of  cloud  and  aerosol  layers  may  be  determined  (Nakajima 
and  King,  1990).  The  high  resolution  of  radiance  data  will  help  in  studying  the  spatial 
structure  of  cloud  systems. 

1.3  Research  Objective 

The  SSP  was  constructed  to  measure  the  spectral  and  polarimetric  reflected  sunlight. 
The  objective  of  this  thesis  is  to  assess  the  ability  of  the  SSP  to  measure  the  spectral  sky 
polarization  to  include: 
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1.  Perform  an  initial  calibration  of  the  SSP  to  enable  the  observation  of  sky  light 
polarization. 

2.  Evciluate  the  performance  of  the  SSP  in  measuring  the  sky  light  polarization. 

Part  1  was  broken  into  three  steps;  installation  of  the  polarization  cubes,  charac¬ 
terization  of  the  optics,  and  determining  the  detector  response  to  incident  light.  Part  2 
involves  the  design  of  an  experiment  to  measure  the  sky  polarization  with  the  SSP,  using 
a  polarized  radiative  transfer  model  to  study  the  effects  local  conditions  have  on  skylight 
polarization,  and  comparing  the  SSP  measured  sky  light  polarization  with  that  predicted 
by  the  radiative  transfer  model. 

1.4  Thesis  Outline 

Chapter  2  describes  the  design  and  specifications  of  the  SSP  and  its  data  acquisition 
system.  In  Chapter  3  the  optics  installation  and  calibration  procedures  are  discussed.  The 
model  used  to  predict  the  sky  light  intensity  and  polarization  is  described  in  Chapter  4 
as  well  as  the  effects  of  Rayleigh  optical  depth,  surface  albedo,  solar  zenith  angle,  and 
haze  loading.  In  Chapter  5  the  measured  sky  light  polarization  is  compared  to  the  model 
predicted  values. 


Chapter  2 


INSTRUMENT  DESIGN 

The  Scanning  Spectral  Polarimeter  (SSP)  is  a  third  generation  radiometer  which 
utilizes  a  rotating  optical  bandpass  filter  to  measure  the  spectral  region  from  400  nm 
through  4000  nm  with  Half  Bandwidths  (HBW)  of  less  than  17  nm  to  60  nm.  Earlier 
versions  of  this  instrument  are  described  by  Stephens  and  Scott  (1985)  and  Scott  and 
Stephens  (1985).  Six  optical  channels  allow  the  measurement  of  flux,  unpolarized  radiance, 
and  the  four  polarized  radiances;  parallel,  perpendicular,  right  hand  circular,  and  left  hand 
circular.  The  optical  layout  of  the  instrument,  a  brief  discussion  of  the  electronics  and 
control  systems,  and  the  data  rate  and  its  format  are  described  in  this  chapter. 

2.1  Radiometer  Overview 

A  photograph  of  the  instrument  is  shown  in  Figure  2.1  and  the  instrument  layout 
is  shown  in  Figure  2.2.  The  instrument  is  composed  of  3  main  components,  the  motor 
drive  assembly  (A),  a  filter  wheel  and  detector  assembly  (B),  and  the  optics  assembly  (C). 
The  optics  are  held  in  a  vacuum  to  reduce  the  hazard  of  thermal  shock  and  condensation 
problems.  The  modular  feature  of  the  instrument  design  offers  a  number  of  advantages, 
including  the  ability  to  use  different  optical  assemblies  (say  with  or  without  polarization 
optics)  or  different  filter  arrangements.  A  more  detailed  discussion  of  each  of  these  main 
components  is  now  presented. 

2.2  The  Optics  Assembly 

An  expanded  view  of  the  optical  configuration  of  the  SSP  is  shown  in  Figure  2.3. 
Radiation  enters  the  instrument  through  one  of  6  windows  which  are  more  clearly  shown 
in  Figure  2.1,  passes  through  the  optical  head  where  it  is  focused  onto  the  Circular  Variable 


Figure  2.1:  Scanning  spectral  polarimeter  instrument  head.  The  six  windows  through 
which  radiation  enters  the  instrument  head  are  shown  on  the  right.  The  valve  on  the  top 
of  the  instrument  is  for  evacuating  the  vacuum  chamber.  Cables  attach  to  the  back  of  the 
instrument  shown  on  the  left.  Overall  dimensions:  6  inch  diameter  and  10.5  inches  long. 

Filter  (CVF)  wheel,  then  through  the  field  stop,  and  onto  the  detector  assembly.  The  field 
stop  is  designed  to  underfill  the  detectors.  The  polarization  and  focusing  optics  are  held  in 
an  optical  tower,  which  allows  each  of  the  optics  to  be  individually  positioned.  Radiation 
of  different  wavelengths  will  focus  at  different  places  since  the  index  of  refraction  varies 
with  wavelength.  Achromatic  focusing  lenses  are  used  to  reduce  this  effect.  Gian- Taylor 
polarization  cubes  axe  used  instead  of  the  more  common  Glan-Thompson.  Glan-Taylors 
have  the  same  optical  configuration  as  the  Glan-Thompson,  except  Glan-Taylors  are  air 
spaced  instead  of  being  cemented  together.  This  reduces  the  thermal  shock  hazards. 

The  optics  of  the  two  circular  polarization  channels  (5  and  6)  axe  characterized  by 
a  fused  silica  window  through  which  radiation  enters  the  instrument.  Radiation  then 
travels  through  a  Fresnel  Rhomb  and  a  Glan-Taylor  polarization  cube.  The  Fresnel  Rhomb 
converts  circular  polarized  radiation  into  linear  polarized  radiation  oriented  at  45°  to  the 
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Figure  2.2:  SSP  components:  the  motor  drive  assembly  (A),  a  filter  wheel  and  detector 
assembly  (B),  and  the  optics  assembly  (C). 

plane  of  incidence.  The  Glan-Taylor  polarization  cube  transmits  only  the  radiation  that  is 
parallel  to  its  optical  axis  and  this  cube  is  oriented  45°  to  the  left  of  the  Fresnel  Rhomb  for 
detection  of  left  hand  circular  polarized  radiation  and  45°  to  the  right  for  detection  of  right 
hand  circular  polarized  radiation.  Radiation  then  travels  through  an  achromatic  focusing 
lens  (made  of  BK7  optical  glciss),  a  second  fused  silica  window  as  it  exits  the  vacuum 
chamber,  is  focused  onto  the  CVF,  passes  through  the  aperture  and  onto  a  detector. 

The  optics  for  the  linear  polarization  channels  (2  and  4)  are  defined  by  entering  a  fused 
silica  window  and  a  Glan-Taylor  polarization  cube  so  aligned  to  allow  only  perpendicular 
or  parallel  polarized  light  to  pass.  The  latter  is  achieved  by  rotating  the  Glan-Taylor 
polarizing  cube  90°.  Radiation  then  travels  through  a  BK7  achromatic  focusing  lens, 
through  a  second  fused  silica  window  to  exit  the  vacuum  chamber,  focuses  on  the  CVF, 
then  through  the  aperture  and  onto  the  detector. 
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Figure  2.3:  SSP  optics. 

The  radiance  channel  (3)  has  a  multi-spectral  Zinc  Sulfide  (ZnS)  window  as  the  en¬ 
trance  port  to  the  instrument,  radiation  then  passes  through  a  ZnS  focusing  lens  at  the 
bottom  of  the  vacuum  chamber  and  exits  through  a  second  ZnS  window,  focuses  on  the 
CVF,  then  through  the  aperture  and  onto  a  detector.  Multi-spectral  ZnS  (also  known  as 
Cleartran)  windows  and  lens  are  used  for  the  path  of  this  channel  due  to  its  fiat  transmis¬ 
sion  curve  (approximately  70%  transmission)  from  400  nm  -  9000  nm.  The  hemispheric 
fiux  channel  (1)  is  similar  to  the  radiance  channel  except  that  radiation  enters  through  a 
flashed  opal  diffuser  window,  the  entrance  and  exit  windows  are  made  of  fused  silica,  and 
the  achromatic  focusing  lens  is  made  of  BK7. 


The  radiance  and  polarization  channels  (2  -  6)  all  have  narrow  fields  of  view  with  a 
viewing  half  angle  of  approximately  20mRad.  The  diffuser  in  channel  1  provides  a  full 
hemispheric  view. 
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2.3  The  Filter  and  Detector  Assembly 

The  CVF  is  made  of  four  90®  section  bandpass  filters  with  a  total  bandpass  of  400 
nm  -  4000  nm.  Table  2.1  shows  the  nominal  spectral  region  covered  by  each  section  and 
the  measured  spectral  region  for  the  filter  sections  used  in  this  experiment.  The  CVFs 
were  commerci2Llly  available.  Variance  between  the  nominal  and  actual  wavelengths  is 
as  much  at  15%,  and  variance  between  individual  filters  with  the  same  nominal  spectral 
range  can  be  as  much  as  5%  of  the  central  wavelength. 


Figure  2.4:  CVF  wheel.  The  optical  encoder  holes  and  detector  positions  are  shown.  The 
detectors  are  stationary,  mounted  on  a  computer  board  below  the  CVF  wheel.  The  CVF 
wheel  rotates  counter  clock  wise. 


The  CVF  sections  are  mounted  on  an  aluminum  wheel  with  a  0.040  inch  space  between 
each  segment.  The  space  provides  a  clear  aperture  for  the  broadband  measurements.  Index 
holes  located  on  the  outside  rim  of  the  aluminum  wheel  specify  points  on  the  filter  wheel 
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Table  2.1:  Nominal  and  spectral  regions  covered  by  each  section  of  the  circular  variable 
filter  (CVF). 


Section 

Nominal 

Spectral 

Range 

Measured 

Spectral 

Range 

1 

400  -  700  nm 

373  -  753  nm 

2 

700  -  1235  nm 

606  -  1384  nm 

3 

1235  -  2225  nm 

1135  -  2235  nm 

4 

2225  -  4000  nm 

2143  -  4091  nm 

where  radiation  data  are  read  by  the  data  acquisition  system,  this  provides  a  custom 
optical  encoder.  A  diagram  of  the  filter  sections  and  encoder  hole  locations  is  shown 
in  Figure  2.4.  Encoder  holes  were  originally  placed  to  at  5.3°  intervals  along  the  visible 
section  (400  -  700  nm)  to  eliminate  overlap  between  adjacent  measurements,  however,  since 
the  detectors  axe  evenly  spaced  around  the  sensor  head,  this  sampling  pattern  wais  only 
valid  for  detector  6.  The  CVF  sections  used  provided  spectral  overlap  between  consecutive 
sections.  The  central  wavelengths  are  listed  in  Appendix  A  for  all  holes  on  the  CVF  as 
a  function  of  index  hole  number  and  position  around  the  filter  wheel.  Table  2.2  lists  the 
central  wavelengths  from  390  -  1100  nm  in  ascending  order.  Bandwidths  for  each  channel 
have  not  been  measured,  but  are  estimated  at  less  than  or  equal  to  4%  of  the  central 
wavelength  for  sections  1  and  2,  and  less  than  1.5%  of  the  central  wavelength  for  sections 
3  and  4. 

2.3.1  Detectors 

The  detectors,  shown  in  Figure  2.5,  are  hybrid  two  color  detectors  with  built  in 
pre~amplifiers.  A  Silicon  (Si)  element  is  mounted  over  a  Lead  Selenide  (PbSe)  element. 
The  Si  element  is  used  in  the  photovoltaic  mode  and  the  PbSe  element  is  used  in  the 
photoconductive  mode.  The  detector  windows  are  also  made  of  multispectral  zinc  sulfide 
(ZnS).  The  Si  detector  has  a  spectral  range  of  400  nm  to  1100  nm  with  an  8  nano¬ 
second  response  time.  The  Si  element  has  a  50%  transmission  at  wavelengths  greater 
than  1100  nm.  This  allows  radiation  to  pass  through  the  element  onto  the  PbSe  detector 
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Table  2.2:  Central  wavelengths  (nm)  in  ascending  order  for  channels  1-4.  *  -  Measure¬ 
ments  where  channels  2  -  4  all  have  the  same  centred  wavelength. 


Value 

Ch  1 

Ch  2 

Ch3 

Ch4 

Ch  4 

1 

394 

394 

394 

24 

687 

666.5 

666.5 

714 

2 

407 

407 

407 

411 

25 

696 

680 

680 

*719.5 

3 

420 

420 

420 

432.5 

26 

712.5 

687 

687 

*733 

4 

433 

433 

433 

453.5 

27 

718 

693 

693 

758 

5 

445 

445 

445 

475 

28 

737.5 

706 

706 

803.5 

6 

458 

458 

458 

*496 

29 

761.5 

712.5 

712.5 

848.5 

7 

470.5 

470.5 

470.5 

30 

788.5 

719.5 

719.5 

896 

8 

483 

483 

483 

*521.5 

31 

814.5 

733 

733 

943 

9 

496 

496 

496 

*535 

32 

842 

737.5 

737.5 

986 

10 

509 

509 

509 

*548 

33 

867 

761.5 

761.5 

1037 

11 

521.5 

521.5 

521.5 

*561 

34 

892 

788.5 

788.5 

1082 

12 

535 

535 

535 

*575 

35 

918 

814.5 

814.5 

13 

548 

548 

548 

*588 

36 

944 

842 

842 

14 

561 

561 

561 

37 

971 

867 

867 

15 

575 

575 

575 

38 

993 

892 

912 

16 

588 

588 

588 

*627 

39 

1019 

919 

958 

17 

601 

601 

601 

40 

1041 

945 

1008 

18 

614 

614 

614 

*653.5 

41 

1077 

970 

1052 

19 

627 

627 

627 

*666.5 

42 

996 

1098 

20 

639 

639 

639 

43 

1025 

21 

650 

640 

640 

44 

1050 

22 

664 

653.5 

653.5 

*693 

45 

1077 

23 

673 

664 

664 

Hiljy 

below.  The  PbSe  detector  has  a  spectral  range  of  1000  nm  to  5200  nm  with  a  response 
time  of  1-3  micro-seconds.  The  upper  wavelength  limit  of  the  PbSe  element  is  temperature 
dependent. 

2.3.2  Sensor  Head  Electronics 

Figure  2.6  shows  a  block  diagram  of  the  SSP  sensor  head  electronics.  The  circular 
variable  filter  is  driven  by  a  computer  controlled  stepper  motor.  Rotation  rate  can  be 
varied  from  0.1  revolution/second  to  30  revolutions /second.  Two  optical  switches  are 
used,  one  to  determine  the  start  of  a  new  scan  and  the  second  generates  data  valid  signals 
from  the  index  holes  on  the  rim  of  the  CVF.  Signals  from  each  of  the  6-2  color  detectors 
are  amplified  so  the  resulting  analog  signal  sent  to  the  data  acquisition  system  varies 
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Figure  2.5:  SSP  two  color  detector.  The  Silicon  element  is  the  black  square  at  the  center  of 
the  detector.  The  Lead  Selenide  detector  has  two  elements,  a  blind  element  shown  as  the 
white  square  below  the  Silicon  element  and  an  active  element  under  the  Silicon  element. 
Amplifiers  are  mounted  on  the  grid  below  the  detectors.  The  entire  setup  fits  in  a  TO-8 
housing  with  overall  dimensions  of  .6  inch  diameter  and  .3  inch  tall. 

from  0  volts/32500  counts  (no  signal)  to  -10  volts/0  counts  (max  signal).  Temperature  is 
monitored  at  four  locations:  the  polarization  and  focusing  optics,  the  CVF  air  chamber, 
the  detector  block,  and  the  motor. 

2.4  Data  Acquisition 

The  SSP  is  designed  to  use  two  different  data  acquisition  systems  (DAS)  to  collect 
information  from  the  sensor  head.  Both  systems  use  a  micro  controller  that  runs  all  timing 
and  data  handling  issues.  Analog  signals  are  received  from  the  sensor  head  on  16  different 
channels.  Descriptions  of  the  16  data  channels  are  listed  in  Table  2.3.  Each  of  the  16 
data  channels  has  its  own  sampling,  16  bit  analog  to  digital  converter.  After  the  A-D 
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Figure  2.6:  Sensor  head  electronics. 


conversion,  digital  data  is  put  into  two  parallel  first  in  first  out  (FIFO)  outputs.  Two 
FIFO  ports  are  provided  to  allow  two  computer  systems  to  read  the  same  data  from  the 
SSP  data  acquisition  system  without  rewiring.  A  real  time  clock  is  used  to  flag  each  data 
scan  with  the  acquisition  time.  The  data  pattern  from  the  FIFO  is  shown  in  Table  2.4 

The  first  DAS  is  rack  mounted  and  was  designed  for  manned  airborne  platforms, 
calibration,  and  ground  testing.  A  block  diagram  of  this  system  is  shown  in  Figure  2.7. 
This  system  has  full  diagnostic  outputs,  real  time  graphical  display,  and  large  data  storage 
capacity. 

The  second  DAS,  a  small  footprint  system,  was  designed  for  unmanned  aerospace 
vehicles  (UAV)  or  experiments  where  unattended  operation  is  required.  This  system 
provides  no  graphical  output,  storage  capacity,  or  input  voltage.  The  user  must  provide 
the  input  voltage,  download  data  from  the  FIFO  ports,  and  transmit  or  store  data  from  the 
sensor  head.  A  block  diagram  of  this  system  is  shown  in  Figure  2.8.  This  small  footprint 
version  does  provide  diagnostic  output  for  an  external  computer  to  aid  in  debugging. 
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Table  2.3:  Data  channel  definitions. 


Data 

Det. 

Wavelength 

Optical 

Ch. 

Det. 

Type 

Region 

Properties 

1 

1 

Si 

400-1100  nm 

Hemispherical 

2 

2 

Si 

400-1100  nm 

Parallel  Pol. 

3 

3 

Si 

400-1100  nm 

Radiance 

4 

4 

Si 

400-1100  nm 

Perpendicular  Pol, 

5 

5 

Si 

400-1100  nm 

Left  Circular  Pol. 

6 

6 

Si 

400-1100  nm 

Right  Circular  Pol. 

7 

1 

PbSe 

1000-4000  nm 

Hemispherical, 

8 

2 

PbSe 

1000-4000  nm 

Parallel  Pol. 

9 

3 

PbSe 

1000-4000  nm 

Radiance 

10 

4 

PbSe 

1000-4000  nm 

Perpendicular  Pol. 

11 

5 

PbSe 

1000-4000  nm 

Left  Circular  Pol. 

12 

6 

PbSe 

1000-4000  nm 

Right  Circular  Pol. 

13 

Detector  Block  Temperature 

14 

Motor  Temperature 

15 
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16 

Polarizing  and  Focusing  Optics  Temperature 
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Figure  2.7:  Rack  mounted  data  acquisition  system. 


Table  2.4;  FIFO  byte  pattern. 


Scan  Header 

First  of  107  Data  Sets 

Byte 

Description 

Byte 

Description 

1 

Clock  (0.01  Seconds) 

27 

Ch  1  (Low  Byte) 

2 

Clock  (Seconds) 

28 

Ch  2  (Low  Byte) 

3 

Clock  (Minutes) 

29 

Ch  3  (Low  Byte) 

4 

Clock  (Hours) 

30 

Ch  4  (Low  Byte) 

5 

Gain  Ch  1 

31 

Ch  5  (Low  Byte) 

6 

Gain  Ch  2 

32 

Ch  6  (Low  Byte) 

7 

Gain  Ch  3 

33 

Ch  7  (Low  Byte) 

8 

Gain  Ch  4 

34 

Ch  8  (Low  Byte) 

9 

Gain  Ch  5 

35 

Ch  9  (Low  Byte) 

10 

Gain  Ch  6 

36 

Ch  10  (Low  Byte) 

11 

Gain  Ch  7 

37 

Ch  11  (Low  Byte) 

12 

Gain  Ch  8 

38 

Ch  12  (Low  Byte) 

13 

Gain  Ch  9 

39 

Ch  1  (High  Byte) 

14 

Gain  Ch  10 

40 

Ch  2  (High  Byte) 

15 

Gain  Ch  11 

41 

Ch  3  (High  Byte) 

16 

Gain  Ch  12 

42 

Ch  4  (High  Byte) 

17 

Motor  Speed  (Low  Byte) 

43 

Ch  5  (High  Byte) 

18 

Motor  Speed  (High  Byte) 

44 

Ch  6  (High  Byte) 

19 

Temp  Ch  13  (Low  Byte) 

45 

Ch  7  (High  Byte) 

20 

Temp  Ch  14  (Low  Byte) 

46 

Ch  8  (High  Byte) 

21 

Temp  Ch  15  (Low  Byte) 

47 

Ch  9  (High  Byte) 

22 

Temp  Ch  16  (Low  Byte) 

48 

Ch  10  (High  Byte) 

23 

Temp  Ch  13  (High  Byte) 

49 

Ch  11  (High  Byte) 

24 

Temp  Ch  14  (High  Byte) 

50 

Ch  12  (High  Byte) 

25 

Temp  Ch  15  (High  Byte) 

26 

Temp  Ch  16  (High  Byte) 

Chapter  3 


INSTRUMENT  CHARACTERIZATION  AND  CALIBRATION 

3.1  Experiment  Description 

For  the  experiments  described  in  this  thesis,  only  four  channels  were  used;  the  ra¬ 
diance,  irradiance,  perpendicular  polarization,  and  parallel  polarizations.  These  channels 
measured  radiation  only  from  400  -  1100  nm  (data  from  the  silicon  detector)  and  the 
broadb2uid  holes  at  each  end  of  the  CVF  section  1  (400  -  700  nm)  were  covered. 

The  first  step  is  to  install  and  align  the  polarization  cubes  in  the  SSP  optical  tower. 
The  SSP  must  then  be  calibrated  using  a  known  light  source.  However,  absolute  calibration 
is  not  crucial  for  the  experiments  described  here,  since  we  are  concerned  with  the  change 
in  sky  polarization  as  a  function  of  view  angle.  Polarization  is  calculated  by  a  ratio  of 
radiance  values,  so  an  accurate  relative  calibration  will  yield  accurate  polarization  values. 
The  pattern  of  sky  polarization  can  then  be  measured  and  compared  to  the  calculated 
pattern  for  a  Rayleigh  atmosphere.  These  steps  will  now  be  discussed. 

3.2  Optics  Alignment 

The  Glan-Taylor  polarization  cube  in  the  perpendicular  channel  must  be  aligned  90° 
from  the  Glan-Taylor  in  the  parallel  channel.  The  equipment  layout  for  this  alignment 
exercise  is  illustrated  in  Figure  3.1.  A  red  unpolarized  lOmW  Helium-Neon  laser  was  used 
as  the  light  source.  The  laser  beam  was  reflected  off  two  movable  mirrors  which  allow  easy 
adjustment  of  the  laser  beam.  The  first  mirror  is  used  to  adjust  the  originating  point  of 
the  laser  beam  as  it  reflects  off  the  second  mirror  and  passes  through  the  optics  tower. 
The  second  mirror  adjusts  the  beam  direction.  The  SSP  optical  tower  was  mounted  on  a 
jig  that  was  built  to  allow  the  optical  tower  to  be  rotated  precisely  180°.  Optical  channels 
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Figure  3.1:  Equipment  setup  for  Glan-Taylor  polarization  cube  alignment. 

2  and  4  are  used  for  the  perpendicular  and  parallel  polarization  so  the  Glan-Taylor  for 
channel  4  is  in  the  same  location  as  the  Glan-Taylor  for  channel  2  after  the  180®  rotation. 
A  silicon  photo  detector  was  mounted  behind  the  optical  tower  to  measure  the  radiation 
transmitted  through  the  calibration  optics.  Lens  centering  guides  were  built  from  a  piece 
of  aluminum  turned  down  to  the  same  diameter  as  the  lens  with  a  1/8  inch  hole  in  the 
center. 

To  install  the  optics,  all  equipment  was  aligned  with  the  laser.  A  reference  polarization 
cube  was  installed  in  lens  holder  1  that  allowed  only  perpendicular  polarized  radiation  to 
pgtss.  The  polarization  cube  for  channel  2  was  installed  in  the  optics  tower  90°  to  this 
reference  cube.  The  SSP  Optics  tower  was  rotated  180°.  A  second  reference  cube  was 
installed  in  lens  holder  2,  90°  from  the  first.  The  polarization  cube  for  channel  4  was 
installed  90°  from  the  second  reference  cube.  Final  alignment  of  the  polarization  cubes 
was  90°  ±  .183°.  The  exact  procedure  used  for  this  exercise  is  described  further  in 
Appendix  B. 
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3.3  Calibration 

To  enable  data  from  the  SSP  to  be  accurately  related  to  atmospheric  spectrum,  the 
instrument  must  be  calibrated  to  a  standardized  source.  For  this  exercise  a  1000  watt 
Quartz  Halogen  DC  light  bulb  was  used.  The  bulb  was  calibrated  by  Eppley  Lab  over  the 
range  of  250  nm  to  2400  nm,  traceable  to  NIST  standards.  The  calibration  curve  for  the 
lamp  is  shown  in  Figure  3.2.  Of  note  are  the  low  irradiance  in  the  spectral  region  below 
500  nm  and  the  irradiance  peak  at  900  nm.  These  characteristics  appear  in  the  following 
calibration  experiments. 


Three  different  calibrations  of  the  SSP  were  performed;  determining  the  field  of  view 
(for  channels  2,  3,  and  4),  measuring  the  angular  response  of  transmission  for  the  flashed 
optical  diffuser  in  channel  1,  and  determining  the  detector  response  to  incident  spectral 
intensity.  Procedures  and  results  of  these  calibrations  will  now  be  discussed. 


Figure  3.3:  FOV  measurement:  Equipment  layout. 

3,3.1  Field  of  View 

The  equipment  setup  used  to  determine  the  field  of  view  for  channels  2  -  4  is  illustrated 
in  Figure  3.3.  This  experiment  was  conducted  in  a  hallway  where  the  distance  between 
the  SSP  and  the  light  source  was  20  feet.  The  SSP  Wcts  set  on  a  leveling  bench  which 
consisted  of  a  thick  aluminum  plate  with  a  height  adjusting  screw  in  each  corner.  The 
light  source  is  on  scissors  jack  that  allows  the  height  of  the  bulb  to  be  estsily  adjusted.  The 
light  source  and  jack  are  on  a  rolling  bench  with  the  wheels  locked  so  the  bench  will  only 
roll  in  a  straight  line.  All  equipment  is  leveled  and  aligned  down  the  center  of  the  hallway 
using  the  Helium-Neon  laser.  A  cardboard  screen  with  a  5  by  7  inch  hole  was  attached  to 
the  rolling  table.  This  screen  moves  with  the  light  source  and  stops  most  of  the  reflection 
from  the  walls.  To  reduce  reflections  further,  two  black  felt  curtains  were  hung  across 
the  hallway  with  a  hole  cut  into  each  curtain  to  allow  the  light  from  the  source  to  reach 
the  SSP.  A  black  felt  curtain  is  also  hung  behind  the  SSP  to  reduce  the  reflection  from 
the  back  wall.  The  source  light  is  moved  perpendicular  to  the  experiment  centerline  as 
shown  in  Figure  3.3.  The  distance  from  the  center  line  and  the  signal  intensity  received 
at  the  SSP  are  recorded.  The  exact  procedure  used  for  this  exercise  is  described  further 
in  Appendix  B. 

Results  of  measured  intensity  versus  light  source  position  for  the  radiance  channel  (3), 
are  shown  in  Figures  3.4.  This  graph  shows  the  relative  power  received  (digital  counts) 
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Figure  3.4:  FOV  measurement:  Relative  intensity  versus  position  and  wavelength  for  the 
radiance  channel  (3).  Distance  between  the  SSP  and  light  source  is  20  ft.  Data  shown  is 
from  the  silicon  detector.  The  jump  in  received  intensity  at  700  nm  is  due  to  CVF  section 
2  having  a  higher  transmission  than  section  1. 


versus  the  light  source  position  and  the  measured  wavelength.  As  the  position  of  the 
source  is  changed  we  see  a  peak  in  the  received  power  at  the  experiment  centerline  with 
the  power  falling  off  rapidly  as  the  light  source  is  moved  off  center.  Looking  at  increasing 
wavelengths,  we  see  an  increase  in  the  power  received  from  400  nm  through  700  nm.  At 
700  nm  there  is  a  large  jump  in  the  power  received,  followed  by  a  steady  increase  to  the 
received  power  maximum  at  900  nm.  Beyond  1000  nm  the  power  received  falls  off  rapidly. 

The  jump  in  received  power  at  700  nm  is  a  factor  of  two  effects:  the  second  CVF 
section  is  physically  thinner  and  therefore  has  a  higher  transmission,  it  also  has  a  steeper 
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wavelength  gradient  and  each  reading  will  have  a  wider  bandwidth.  The  rapid  decrease 
in  received  power  beyond  1000  nm  is  a  combined  factor  of  the  lamp  emitting  less  in  this 
spectral  region  and  approaching  the  effective  range  of  the  silicon  detector.  Figure  3.5 
shows  the  corresponding  graph  for  the  parallel  polarization  channel  (2).  Features  of  the 
parallel  polarization  are  similar  to  the  radiance  channel,  except  the  power  received  is 
approximately  half  that  of  the  radiance  channel. 


Figure  3.5:  FOV  measurement:  Relative  intensity  versus  position  and  wavelength  for  the 
parallel  polarization  channel  (2).  Distance  between  the  SSP  and  light  source  is  20  ft.  Data 
shown  is  from  the  silicon  detector.  The  jump  in  received  intensity  at  700  nm  is  due  to 
CVF  section  2  having  a  higher  transmission  than  section  1. 


22 


Figure  3.6:  FOV  measurement:  Normalized  intensity  versus  position  and  wavelength  for 
the  radiance  channel.  Distance  between  the  SSP  and  light  source  is  20  ft.  Data  shown  is 
from  the  silicon  detector. 


The  field  of  view  is  determined  by  finding  the  two  positions  where  the  received  power 
is  ^  the  maximum.  The  width  between  these  two  positions  and  the  distance  between  the 
light  source  and  the  SSP  are  then  used  to  calculate  the  viewing  angle.  To  automate  the 
process,  the  data  were  first  normalized  as  shown  in  Figure  3.6  for  the  radiance  channel. 
Note  how  the  constant  distant  lines  roll  over  the  top  of  the  curve  at  the  near  infra  red 
wavelengths.  This  is  due  to  the  longer  wavelength  light  being  ’bent’  ats  is  passes  through 
the  optics.  In  effect  the  SSP  looks  in  a  slightly  different  direction  at  the  longer  wavelengths. 

The  data  for  each  wavelength  are  now  fit  to  a  gaussian  curve  of  the  form: 


Normalized  Power 


23 


-40  -30  -20  -10  0  10  20  30  -40  -30  -20  -10  0  10  20  30  40  ’ 


Source  Position  (cm) 

Figure  3.7:  FOV  measurement:  Gaussian  curve  fit  to  the  normalized  data.  Data  for  521 
nm  is  shown  in  (A)  and  912  nm  in  (B).  Distance  between  the  SSP  and  light  source  is  20 
ft. 

k 

ReceivedPower  =  —  exp[— .5{(d  —  (3.1) 

s 

Where  d  =  distance  from  the  experiment  centerline,  m  is  the  distribution  mean,  s  is  the 
standard  deviation,  and  k  is  the  amplitude  factor.  Figure  3.7  shows  how  the  actual  data 
fit  the  gaussian  curve  for  521  nm  and  912  nm.  Once  the  parameters  k^  s,  and  in  are  known 
the  equation  can  be  solved  for  the  two  locations  where  received  power  equals  50%  of  the 
maximum,  or  more  appropriately,  the  width  of  the  50%  viewing  area 

ViewWidth  =  y(2m)2  -  4(m2  -|-  s‘^{2ln{^)))  (3.2) 

where  p  is  50%.  Once  the  the  viewing  width  is  known,  it  is  trivial  to  calculate  the  half 
angle  of  view.  The  resultant  half  angle  fields  of  view  are  shown  in  Figure  3.8. 

3,3.2  Cosine  Response  of  Diffuser 

The  object  of  this  exercise  is  to  measure  how  the  transmission  of  the  flashed  opal 
diffuser  in  channel  1  varies  as  a  function  of  the  incident  angle  of  incoming  radiation. 
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Wavelength  (nm) 

Figure  3.8:  Field  of  view  versus  wavelength  for  the  radiance,  parallel  polarization,  and 
perpendicular  polarization  channels. 

Flashed  opal  is  the  only  commercially  available  material  that  allows  the  transmission  of 
radiation  from  any  angle.  The  diffuser  is  made  by  bonding  a  thin  layer  of  diffusing  glass 
to  a  supporting  clear  glass  substrate.  The  diffusing  glass  causes  multiple  scattering  of 
light.  By  nature,  the  shorter  wavelengths  scatter  more,  which  also  means  transmission 
losses  are  greater.  Flashed  opal  diffusers  work  well  through  visible  wavelengths  and  the 
transmission  curve  is  close  to  that  shown  in  Figure  3.9,  however,  at  longer  wavelengths  a 
directed  component  is  superimposed  at  the  forward  direction.  Part  of  this  exercise  is  to 
define  the  extent  of  this  forward  transmission  lobe. 

A  schematic  of  the  diffuser  response  setup  is  shown  in  Figure  3.10.  Setup  is  similar  to 
that  for  the  narrow  field  of  view  exercise,  however,  the  distance  between  the  SSP  and  the 
light  source  is  only  53.3  cm.  The  lamp  was  calibrated  by  Eppley  Labs  at  a  distance  of  50 
cm,  they  recommend  use  of  the  lamp  only  for  distances  between  the  lamp  and  detector  in 
excess  of  45  cm.  The  SSP  was  set  as  close  as  it  could  be  to  the  lamp  and  still  have  room 
for  the  screens  to  reduce  reflections.  Once  the  SSP  and  light  source  were  aligned,  the  SSP 
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Double  Ground 


Figure  3.9:  Transmission  of  an  ideal  diffuser,  a  flashed  opal  diffuser,  and  double  ground 
quartz. 


was  rotated  and  the  received  intensity  was  measured  as  the  angle  between  the  SSP  and 
light  source  was  varied.  The  procedure  used  is  described  in  detail  in  Appendix  B. 


Figure  3.10:  Diffuser  transmission:  Equipment  setup  for  measuring  the  angle  dependence 
of  transmission  for  the  flashed  opal  diffuser. 


Results  of  the  measured  intensity  versus  the  incident  angle  are  shown  in  Figure  3.11. 
Data  were  not  recorded  when  angles  exceeded  80°  because  the  front  of  the  SSP  vacuum 
chamber  blocked  all  light  from  reaching  the  diffuser  lens.  In  general,  power  increases  as  the 
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Figure  3.11:  Diffuser  Transmission:  Relative  power  versus  wavelength  and  angle  between 
SSP  and  Source.  Distance  between  the  SSP  and  light  source  is  53.3  cm. 

source  angle  is  decreased  to  0.  From  400  nm  through  800  nm  there  is  a  gradual  increase, 
although  the  received  power  is  so  low  it  is  difficult  to  see.  From  800  nm  through  1100  nm 
there  is  an  extreme  increase  in  power  received  at  angles  less  than  3^,  this  is  the  strong 
forward  transmission  lobe  mentioned  previously.  The  maximum  power  received  is  shifted 
slightly  to  1000  nm  instead  of  900  nm.  Beyond  1100  nm  the  signal  falls  off  sharply  as  this 
is  beyond  the  range  of  the  silicon  detector. 

To  further  explore  how  the  flashed  opal  diffuser  transmits  radiation,  the  data  were 
normalized  so  the  maximum  signal  received  is  equal  to  unity  as  shown  in  Figure  3.12. 
Three  distinct  areas  of  diffuser  transmission  are  evident.  From  500  nm  through  750  nm 


Figure  3.12:  Diffuser  transmission:  Normalized  power  versus  wavelength  and  angle  be¬ 
tween  SSP  and  source.  Distance  between  the  SSP  and  light  source  is  53.3  cm. 

the  transmission  curve  is  roughly  a  cosine  curve.  However,  between  500  and  600  nm  the 
received  signal  is  weak  enough  that  instrument  noise  distorts  the  curve.  Below  500  nm, 
the  data  were  to  noisy  too  define  a  curve.  From  850  nm  through  1100  nm  the  transmission 
curve  has  a  very  strong  peak  between  0°  and  2°.  At  1000  nm  the  power  received  drops 
to  less  than  85%  of  the  maximum  by  3^^.  Between  750  nm  and  850  nm  there  is  a  rapid 
transition  from  a  cosine  transmission  curve  to  a  curve  with  the  strong  forward  transmission 


28 


Figure  3.13:  Diffuser  transmission:  Normalized  power  versus  wavelength  and  angle  be¬ 
tween  SSP  and  source,  forward  transmission  lobe  removed.  Distance  between  the  SSP 
and  light  source  is  53.3  cm.  Forward  transmission  lobe  is  removed  by  eliminating  received 
power  data  for  2°  to  0°. 

The  forward  transmission  lobe  wa^  removed  by  normalizing  the  data  between  3°  and 
80®  to  unity  as  shown  in  Figure  3.13.  Even  at  the  longer  wavelengths,  the  transmission 
curve  approximates  a  cosine  curve  once  the  forward  transmission  lobe  is  removed. 

To  allow  the  calculation  of  diffuse  light,  a  formula  for  the  transmission  curve  is  needed. 
For  this  purpose  a  curve  of  the  form: 


ReceivedPower{9)  ^  r 

&cei.erP<,».r(9  =  0)  =  <*«  +  ‘«  +  «  + 


(3.3) 
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Figure  3.14:  Diffuser  transmission:  Curve  fit  for  the  diffuse  channel  transmission  function. 
Data  for  521  nm  is  shown  in  (A)  and  912  nm  in  (B).  Distance  between  the  SSP  and  light 
source  is  53.3  cm. 

where  the  parameters  a,  6,  c,  d,  and  e  are  all  adjusted  to  fit  the  data.  The  first  three  terms 
of  the  right  hand  side  will  approximate  the  cosine  transmission  curve.  The  exponential 
term  accounts  for  the  forward  transmission  lobe.  The  parameter  e  is  adjusted  to  cut  off  the 
effects  of  the  forward  transmission  lobe  at  the  proper  angle.  Figure  3.14  shows  an  example 
curve  fit  for  both  521  nm  and  918  nm.  The  parameters  for  all  the  diffuse  wavelengths  are 
tabulated  in  Appendix  C. 
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3.3.3  Detector  Response 

The  detector  response  was  determined  from  two  tests  -  the  first  applies  to  the  diffuse 
channel  (1)  and  the  second  to  the  three  narrow  field  of  view  channels  (2  -  4)  channels.  For 
the  diffuse  channel  the  SSP  looks  directly  at  the  calibration  light  and  measurements  are 
taken  to  relate  SSP  measured  voltage  counts  to  the  incident  radiation  flux.  For  the  three 
narrow  filed  of  view  channels  the  light  source  is  reflected  from  a  lambertian  reflectance 
plate,  the  SSP  looks  at  the  lambertian  plate.  Again,  a  relationship  between  voltage  counts 
and  incident  radiation  is  determined. 

Data  processing  and  measurement  errors  have  not  yet  been  addressed.  SSP  raw  data 
consist  of  a  stream  of  8  bit  values,  described  in  Table  2.4.  One  scan  of  data  is  normally 
2594  bytes  long,  however,  for  this  experiment  the  data  stream  were  shortened.  Since  data 
from  only  4  of  the  12  detectors  are  used,  the  byte  pattern  was  shortened  to  864  bytes. 
This  decreased  the  amount  of  data  going  to  the  data  acquisition  computer  hard  drive,  and 
fewer  data  write  errors  are  encountered.  Under  normal  conditions,  the  CPU  data  write 
process  lags  the  SSP  data  acquisition  and  approximately  1  scan  of  20  will  not  be  written 
to  disk  correctly.  With  the  shortened  byte  pattern,  the  SSP  normally  drops  one  scan  of 
the  first  5  and  then  accurately  records  all  data.  The  analog  to  digital  converters  also  have 
a  maximum  value  they  will  assign  to  the  high  byte.  Zero  signal  from  the  detectors  is  a 
maximum  digital  count  and  detector  saturation  has  a  zero  digital  count.  If  the  CPU  data 
write  process  gets  behind  the  SSP  acquisition,  the  normal  error  is  to  shift  the  byte  pattern 
one  byte  up  or  down.  To  find  ’bad  data’  the  length  of  each  scan  is  checked  as  it  is  read 
from  disk.  Scans  that  are  too  short  or  too  long  are  considered  bad.  All  high  bytes  are 
also  checked  to  verify  they  do  not  exceed  the  maximum  high  byte  value  due  to  a  data 
shift.  A  scan  with  the  data  shifted  is  also  considered  bad.  Since  there  are  generally  only 
1-2  bad  scans  with  each  data  file,  recovery  of  data  from  the  bad  scans  is  not  critical  and 
the  bad  scans  are  just  thrown  out.  Additionally,  the  first  5  scans  of  each  data  set  are 
thrown  out  to  allow  the  SSP  electronics  to  settle  into  the  normal  operating  mode  after 
the  start  up  power  surge.  After  each  scan  has  been  verified  good,  it  is  added  into  both 
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an  averaging  array  and  an  error  array.  To  calculate  errors  for  the  received  power,  the 
standard  deviation  is  determined  using  the  computational  routine: 


^Vo  = 


f  53^?  "" 

n  —  1 


(3.4) 


Background  noise  causes  a  DC  voltage  offset  to  be  introduced  into  the  data.  This 
offset  varies  for  each  channel  and  each  gain  setting.  Figure  3.15  shows  one  scan  of  raw 
data  .  To  remove  the  voltage  offset,  data  are  normalized  so  readings  from  the  two  covered 
broadband  holes  are  equal  to  zero.  The  data  can  then  be  parsed  and  the  proper  wavelength 
assigned  to  each  reading. 


Diffuse  Channel 


Setup  for  the  detector  response  of  the  diffuse  channel  is  illustrated  in  Figure  3.16. 
The  diffuse  channel  looks  directly  at  the  calibration  source.  Alignment  is  achieved  using 
the  the  laser.  Once  the  SSP  is  aligned  with  the  source  lamp,  received  power  is  measured 
for  every  gain  setting  the  SSP  uses  to  record  data.  The  distance  between  the  SSP  and 
light  source  was  only  51.4  cm.  Measurement  of  the  received  power  is  related  to  the  light 
source  spectral  intensity  by 


kx 


fsourcex  HBWxjdcal/dfTx 
Vc 


(3.5) 


where  k\  is  the  calibration  coefficient,  fsourcex  spectral  intensity  emitted  by  the  cal¬ 

ibration  lamp,  HBWx  is  the  half  bandwidth,  deal  is  the  calibration  distance  for  the  source 
(50  cm),  d  is  the  distance  between  the  SSP  and  source  lamp,  Tx  is  the  path  transmission, 
and  Vc  is  the  measured  voltage  counts.  Due  to  the  short  path  length  and  low  relative 
humidity  during  the  calibration  (less  than  25%)  the  transmission  is  assumed  to  be  1. 


NFOV  Channels 


Setup  for  determining  the  detector  response  of  the  narrow  field  of  view  channels  is 
illustrated  in  Figure  3.17.  The  purpose  of  this  exercise  is  to  relate  the  voltage  counts 
received  from  the  SSP  for  different  channels  and  gain  settings  to  both  each  other  and  the 
incident  radiance.  The  calibration  lamp  illuminates  the  lambertian  reflectance  plate  and 
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Diffuse  and  Radiance  Channels 


Polarization  Channels 


Figure  3.15:  SSP  raw  data.  Gain  setting  on  channel  1  is  25.75.  Gain  setting  for  channels 
2-4  is  13.38.  The  zero  signal  point  is  at  the  top  of  the  graph  and  the  curve  drops  as 
received  power  increases. 


the  SSP  measures  the  reflected  radiance.  The  laser  is  again  user  to  align  the  equipment, 
then  received  power  measurements  are  made  for  each  gain  setting  the  SSP  records  data. 
This  is  the  most  critical  portion  of  the  calibration.  The  setup  procedure  used  is  described 
in  Appendix  B. 

Spectral  intensity  emitted  by  the  source  is  related  to  received  power  by 

fsourcex  H  BW\{^dcal  /  RxT\^_j^Tx^_ggp  COS  6  COS  (j) 


kx  = 


PirVr 


(3.6) 


where  kx  is  the  calibration  coefficient,  fsourcex  ^he  spectral  intensity  emitted  by  the 


calibration  lamp,  HBWx  is  the  half  bandwidth,  deal  is  the  calibration  distance  for  the 
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Figure  3.16:  Equipment  setup  for  determining  the  detector  response  for  the  diffuse  chan- 
nel(l). 

source  (50  cm),  d  is  the  distance  between  the  reflectance  plate  and  source  lamp,  Rx  is 
the  spectral  reflectance  for  the  lambertian  plate,  Tx^_^  is  the  path  transmission  from  the 
source  lamp  to  the  reflectance  plate,  Tx^_g^p  is  the  path  transmission  from  the  reflectance 
plate  to  the  SSP,  cos0  is  the  angle  the  between  the  reflectance  plate  normal  and  the 
incident  source  radiation,  cos(f)  is  the  angle  between  the  SSP  and  the  reflectance  plate 
normal,  P  accounts  for  the  Glan-Taylors  in  the  polarization  channels  and  is  equal  to  1  for 
the  radiance  channel  and  2  for  the  parallel  and  perpendicular  polarization  channels,  and 
Vc  is  the  measured  voltage  counts.  As  long  as  the  lambertian  plate  fills  the  entire  field 
of  view,  the  NFOV  angle  and  the  distance  between  the  SSP  and  lambertian  plate  cancel 
from  the  equation.  Again,  the  path  transmissions  were  assumed  equal  to  1.  Resultant 
calibration  coefficients  enable  measured  voltage  counts  to  be  converted  into  atmospheric 
spectral  intensity. 

To  calculate  the  error  associated  with  the  calibration  coefficients,  the  calibration  equa¬ 
tions  (3.5)  and  (3.6)  must  be  considered.  For  an  absolute  calibration,  every  term  could 
induce  errors.  The  emission  from  the  source  changes  with  power  fluctuations,  the  HBW’s 
vary  with  each  wavelength,  measuring  the  distance  between  the  SSP  and  source  inaccu¬ 
rately  induces  an  error,  path  transmission  changes  with  room  conditions,  and  instrument 
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Figure  3.17:  Equipment  setup  for  determining  the  detector  response  for  the  narrow  field 
of  view  channels  (2-4). 

electronic  noise  causes  the  received  voltage  counts  to  vary.  The  largest  errors  arise  from 
the  instrument  placement  (both  alignment  and  distance  between  SSP  and  light  source) 
and  the  instrument  noise.  Having  the  SSP  out  of  alignment  will  induce  a  large  error  for 
the  diffuse  channel,  especially  at  longer  wavelengths.  A  1  cm  measurement  error  in  the 
distance  between  the  SSP  and  source  causes  a  4%  error  in  the  calibration  equation.  How¬ 
ever,  since  we  are  only  interested  in  a  relative  calibration  between  channels  most  of  these 
errors  cancel.  The  only  error  we  need  to  consider  is  the  instrument  noise,  or  variation  in 
received  voltage  counts  and  the  error  for  the  calibration  coefficient  is 

=  H^).  (3.7) 

The  calibration  coefficients,  A;,  for  the  diffuse  and  radiance  channels  are  shown  in 
Figure  3.18.  The  plot  shows  both  the  averaged  k  values  and  the  associated  error.  Several 
interesting  features  emerge  from  Figure  3.18.  First,  the  calibration  coefficients  for  the 
diffuse  channel  are  about  5  times  that  of  the  radiance  channel.  This  is  due  to  the  low 
transmission  of  the  flcished  opal  diffuser  lens.  Since  less  of  the  incident  radiation  is  able 
to  reach  the  detector,  the  calibration  coefficient  must  be  larger.  Second,  the  error  at  the 
shorter  wavelengths  are  much  larger  than  the  longer  wavelengths.  This  is  a  side  effect  of  the 
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calibration  source.  The  emission  from  the  calibration  source,  shown  in  Figure  3.2,  is  small 
below  500  nm.  At  the  shorter  wavelength  there  is  less  received  power,  and  instrument  noise 
is  a  larger  portion  of  the  measured  signal.  For  the  diffuse  channel,  additional  radiation 
losses  from  the  fiajshed  opal  lens  extend  this  effect  and  the  calibration  coefficients  have 
significant  errors  up  to  700  nm.  Third,  the  calibration  coefficients  increase  rapidly  at  1050 
-  1100  nm.  This  is  the  effective  range  of  the  silicon  detectors,  and  even  though  radiation 
energy  is  reaching  the  detector  element,  they  do  not  measure  it.  Fourth,  the  calibration 
coefficients  for  the  radiance  channel  jump  at  the  end  of  CVF  segment  1  740  nm)  and 

at  the  beginning  of  CVF  section  2  660  nm).  This  is  caused  by  the  cover  put  on  the 

broadband  channel  partially  covering  the  filter  for  these  wavelengths.  The  same  occurs 
for  all  channels,  the  data  for  these  wavelengths  was  removed  from  the  diffuse  channel  plot. 
Fifth,  the  calibration  coefficients  for  CVF  section  2  are  smaller  than  section  1  since  section 
2  has  a  higher  transmission. 


K  Value 


\n 


Figure  3.18:  Calibration  coefficients  for  the  diffuse  (Ch  1)  and  radiance  (Ch  3)  channels. 


Chapter  4 


MODEL 


4.1  Introduction 

A  polarized,  plane  parallel,  monochromatic  radiative  transfer  model  developed  for 
a  previous  activity  was  adapted  for  this  study  (Evans  and  Stephens,  1991).  It  applies 
a  rigorous  solution  of  the  radiative  transfer  equation  for  polarized  light  using  the  four 
Stokes  parameters.  For  clear  sky  radiation,  monochromatic  works  well  since  scattering 
and  absorption  properties  are  uniform  over  the  bandwidths  of  the  SSP  channels.  A  3D 
cartesian  coordinate  system  is  used,  with  the  assumption  of  a  plane  parallel  atmosphere. 
A  point  in  the  medium  is  defined  by  its  2  coordinate  and  the  angular  direction  defined  by 
the  directional  vector  This  vector  is  represented  by  the  angular  pair  (0, 0)  where  0  is 
the  zenith  angle  and  (j)  the  azimuth  angle,  see  Figure  4.1. 


Figure  4.1:  The  general  co-ordinate  system  and  the  specification  of  {0,(j)). 
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This  model  calculates  the  Stokes  parameters  {I,Q,  U,  V)  for  radiation  diffusely  scat¬ 
tered  by  the  atmosphere.  The  complete  polarized  state  of  the  radiation  field  is  described 
by  the  intensity  vector 

l  =  [I,Q,U,V]  (4.1) 

Where  I  is  the  total  intensity,  Q  characterizes  the  degree  of  polarization,  U  the  plane  of 
polarization,  and  V  the  ellipticity  of  the  intensity.  The  Stokes  parameters  have  units  of 
intensity  or  power.  This  allows  the  Stokes  parameters  for  multiple  incoherent  waves  to 
be  added  together  to  obtain  the  total  polarized  state  of  a  radiation  field.  Polarization 
may  range  anywhere  from  unpolarized  (Q  =  t/”  =  F  =  0)  to  totally  polarized  (/^  = 

Q2^i;2^Y2y 

The  Stokes  parameters  may  be  expressed  in  terms  of  the  simple  wave  complex  electric 
field  vector  (Ey^  Eh)  by: 

I=\Evf  +  \EHf 

Q  =  \Ev'^  - 
U  =  2\Ey\\EH\cos6 

V  =  2\Ey\\EH\sin6  (4.2) 

where  6  is  the  phase  difference  between  Ey  and  Eh. 

An  alternative  set  of  Stokes  parameters  that  separates  the  parallel  (J/)  and  perpem 
dicular  (7^)  components  is  defined  by: 

I  =  [7/,/„C/,P]  (4.3) 

where  7/  and  7^  are  the  intensity  components  perpendicular  to  each  other  and  to  the 
direction  of  propagation.  The  SSP  measures  7/  and  Ir  directly. 
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4.2  The  polarized  radiative  transfer  equation 

The  radiative  transfer  equation  that  describes  the  transfer  of  I  through  the  plane 
parallel  atmosphere  can  be  written  as 

M— =  -I(r,  ^  M(m,  <PWd<f>'  +  <r(r,  fi,  4>)  (4.4) 

where,  I  is  the  four  component  Stokes  vector,  M  is  the  four-by-four  scattering  matrix 
which  depends  on  the  Stokes  representation  (4.1),  a  a  vector  for  radiation  sources,  Q  is 
the  single  scatter  albedo,  r  is  the  optical  depth,  /i  is  the  cosine  of  the  zenith  angle,  and  (j) 
is  the  azimuth  angle.  The  detailed  structure  of  M  will  be  discussed  latter.  The  coordinate 
system  used  is  r  increases  downward  and  fi  is  positive  for  downward  directions. 

The  radiation  field  is  separated  into  collimated  and  diffuse  components.  The  above 
radiative  transfer  equation  involves  only  the  diffuse  component.  The  collimated  direct 
solar  beam  is  attenuated  according  to  Beer’s  Law  and  the  radiation  scattered 

from  the  collimated  solar  beam  becomes  a  ’’pseudo-source”. 

The  radiation  sources  can  be  separated  into  two  terms, 

a{T,  fi,  <f>)  =  ci>;  -1^0, 4>oW  +  (1  -  ‘^)B(r)  (4.5) 

The  first  term  on  the  right  hand  side  of  (4.4)  describes  the  source  of  diffuse  I  due  to  single 
scatter  from  the  direct  solar  beam.  Solar  radiation  is  unpolarized  and  of  the  form 

F  =  [F,/mo,  0,0,0]  (4.6) 

where  Fq  is  the  direct  solar  flux  at  the  top  of  the  atmosphere  and  is  taken  to  be  incident 
along  the  direction 

The  second  term  on  the  right  hand  side  of  (4.4)  accounts  for  thermal  emission  of  the 
medium.  This  emission  function  is  also  unpolarized  and  of  the  form 


B  =  [B{T), 0,0,0] 


(4.7) 
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and  B{T)  is  the  Planck  black  body  function 

1.1911  X  10® 


B(T)  = 


L4388^<i0l\  _  n 
AT  ^ 


(4.8) 


X^[exp{- 

where  A  is  the  wavelength  in  microns,  T  is  the  temperature  in  Kelvins,  and  the  units  of 
B  are  Watts/ (meter^  ster  micron). 


4.3  The  phase  matrix 


Single  scattering  calculations  provide  the  scattering  phase  matrix  P(/x,  (f>;  //',  </>')  which 
characterizes  the  single  scatter  of  the  Stokes  parameters  and  is  defined  with  respect  to 
the  scattering  plane,  or  the  plane  defined  by  the  directions  of  incident  and  scattered 
radiation.  Multiple  scattering  calculations  require  that  the  scattering  matrix  be  defined 
with  respect  to  a  fixed  frame  of  reference  rather  than  the  relative  frame  of  reference  which 
defines  P(/i,  /x',  ^').  Thus  a  transformation  is  required  that  accounts  for  this  rotation 

and  which  refers  the  Stokes  parameters  of  the  incident  radiation  to  the  meridional  plane 
containing  the  direction  of  incidence  0')  and  the  Stokes  parameters  of  the  scattered 
radiation  to  the  meridional  plane  containing  the  direction  of  scattering  (/i,<^).  Such  a 
transformation  has  been  discussed  in  detail  by  Chandrasekhar  (1960)  and  is  defined  by 
the  rotation  matrix  L(x) 


L(z) 


/I  0  0  0\ 

0  cos2i  —  sin2z  0 
0  sin2i  cos2i  0 


(4.9) 


\0  0  01/ 


for  the  (I,Q,U,V)  Stokes  representation. 

The  polarization  reference  plane  is  rotated  from  the  incident  meridional  plane  to  the 
scattering  plane,  the  single  scattering  transformation  represented  by  P  is  applied,  and  the 
polarization  is  rotated  to  the  outgoing  meridional  plane  (see  Figure  4.2).  The  result  of 
these  transformations  is  a  phase  matrix  M(/i,/i';^  —  ^')  defined  with  respect  to  a  fixed 
frame  of  reference  where 


M(/z,  /x';  (j) -(}))  =  L(7r  -  i2)P(M,  <t>] 


(4.10) 
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Figure  4.2:  Illustration  of  scattering  angle  and  the  rotation  of  the  polarization  reference 
frame.  The  polarization  of  the  incident  ray  terms  of  the  PiOZ  plane,  is  ro¬ 

tated  into  the  scattering  plane  (P1OP2),  the  scattering  matrix  is  applied,  and  then  the 
polarization  of  the  scattered  ray  (0,  (f))  is  rotated  into  the  P2OZ  plane. 

The  scattering  angle  0  in  Figure  4.2  is  formally  defined  as 

cos©  =  cos  0  cos  0'  +  sin0  +  sin0'cos(</>  —  (/>').  (4.11) 

The  angles  ii  and  22  are  the  angles  that  the  plane  of  scattering  makes  with  the  meridian 
planes  containing  the  directions  of  incidence  and  scattering  respectively.  These  angles  can 


be  determined  from  spherical  trigonometry  by 

sin 21  =  sin^sin(0^  ””0)/ sin©  (4.12) 

sin  22  =  sin  sin((?!>^  —  </>)/ sin©  (4.13) 

cos  21  =  [sin  cos  0  —  sin0cos0^cos(0^  —  (/>)]/ sin©  (4.14) 

cos 2- [sin 0 cos —  sin cos 0 cos (^^  “  (f))]/ sin©  (4.15) 


For  randomly  oriented  particles  with  a  plane  of  symmetry  the  four-by-four  scattering 
phase  matrix  has  only  six  unique  values  (Hovenier,  1969); 


(Pi 

P2 

0 

P2 

P5 

0 

0 

0 

0 

P3 

Pa 

V  0 

0 

-Pa 

Pe) 

P(cos©)  = 


(4.16) 
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The  phase  matrix  for  spheres  has  Pi  =  P5  and  P3  =  Fe- 

After  the  polarization  rotations  have  been  applied,  the  scattering  matrix  has  the  form 


(4.17) 


Pi 

P2  COS  2ii 

— P2sin  2ii 

0  \ 

P2  cos  2t2 

P5COS  2ii  cos  2^2 
— Pasin  2ii  sin  2^2 

— F5sin2ii  cos2t2 
— P3C0S  2ii  sin  2^2 

— F4sin  212 

P2  sin  2^2 

P5COS  2ii  sin  2^2 
+P3sin2ii  cos2z2 

— P5  sin  2ii  sin  2^2 
+P3  cos  2i\  cos  2^2 

P4  cos  2^2 

0 

— P4  sin  2ii 

— P4C0S  2ii 

Pe  / 

4.4  Numerical  Considerations 

4.4.1  Fourier  Expansion  in  Azimuth 

The  first  step  in  formulating  the  solution  to  (4.4)  is  to  represent  the  azimuthal  depen¬ 
dence  of  all  functions  as  a  Fourier  series.  A  real  function  F{(j))  is  expanded  in  the  finite 
series 

M 

F{<t>)  =  5]  cos  ^  sin  4>  F^  =  0  (4.18) 

m— 0 

with  the  sines  and  cosines  forming  and  orthogonal  set  of  functions 

r27r 

/  cos  m(j)  COS  m(j>  = 


rZTT 

/  COS  m(j)  cos  m(j>  = 

Jo 

/  sin  m(f)  sin  m<f)'  =  c^6fmm! 

Jo 

r2ir 

/  COS  mcj)  sin  mcj)^  =  0 


where 


27r  r  =  m  =  0 
^m  —  ^  0  r  =  s,  m  =  0 

TT  m  >  0 


(4.19) 

(4.20) 

(4.21) 


The  integral  over  azimuth  of  a  product  of  functions  transforms  to  a  sum  of  the  product 
of  the  Fourier  coefficients 

f2n  ^ 

/  A{<t>)B{<i>)d{<t>)  =  E  E  (4.22) 

m=0T=c,3 

The  radiative  transfer  equation  is  transformed  to 

dlmr  ^  ^  y  ^  !  /V  /  _ _ _ 


(4.22) 


M  ,  --  Imr  +  .  f  M  mm*rr'  (/^?  M  (/^  )^M  (4.23) 

^  ^  m'=0r*— c, 3 
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m  =  0, r  =  c,  s 

where  the  m  subscript  refers  to  the  azimuthal  Fourier  mode,  and  the  r  subscript  refers  to 
the  cosine  and  sine  terms.  The  Fourier  expansion  is  done  separately  for  each  element  of 
the  Stokes  radiance  vectors  and  the  scattering  matrix. 

4.4.2  Numerical  Quadrature  in  Zenith  Angle 


The  integral  in  zenith  is  approximated  with  a  numerical  quadrature  scheme,  which 
uses  polynomial  interpolation  to  numerically  perform  the  integration.  A  quadrature 
scheme  provides  a  set  of  abscissae  and  weights  for  the  integral  approximation 


f  f{x)dx^j^Wif{Xi) 
1  =  1 


(4.24) 


For  our  application,  the  integral  is  over  the  cosine  of  the  zenith  angle(/2),  or  a  =  —1  to 
6  =  1.  Abscissae  (/Xj)  and  weights  are  determined  for 


ri  N  N 

j=i  i=i 


(4.25) 


where  the  first  and  second  terms  of  the  right  hand  side  represent  the  contribution  to  the 
diffuse  radiation  field  from  the  scatter  over  the  upper  and  lower  hemispheres  respectively. 
There  are  N  angles  in  each  hemisphere  of  solid  angle.  This  approximation  is  sometimes 
referred  to  as  the  ”2n  stream”  approximation. 

The  model  provides  a  choice  of  four  quadrature  schemes:  Gaussian,  double  Gaussian, 
Lobatto,  and  extra  angles.  In  the  Gaussian  scheme,  the  angles  are  provided  by  the  method 
and  are  chosen  so  the  angle  /x  is  not  equal  to  -1,  0,  or  1.  The  Gaussian  scheme  uses  the  same 
angles  in  each  hemisphere.  Gaussian  quadrature  can  exactly  integrate  polynomials  up  to 
degree  2n  —  1.  It  is  exact  for  integration  over  both  hemispheres,  but  not  for  integrations 
over  only  one  hemisphere. 

The  double-Gaussian  is  used  if  accurate  hemispheric  fluxes  are  desired.  This  method 
uses  integration  limits  of  0  to  1  and  creates  a  slight  clustering  of  points  near  one  side  of 
the  interval.  It  works  well  for  low  solar  angles. 

The  Lobatto  quadrature  scheme  is  similar  to  the  Gaussian  method,  only  the  limits 
of  integration  are  included  in  the  abscissae.  This  provides  a  quadrature  angle  in  both  the 
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zenith  and  nadir.  This  method  is  well  suited  for  many  remote  sensing  applications  that 
use  a  zenith  or  nadir  viewing  direction. 

The  extra  angle  mode  uses  the  Gaussian  scheme,  but  allows  the  user  to  select  angles 
to  add  to  the  abscissae  the  method  generates.  The  extra  angles  are  assigned  a  weighting 
factor  of  zero.  This  mode  provides  an  easy  method  of  comparing  the  model  results  to 
tabulated  values. 

Replacing  the  zenith  angle  integral  with  the  quadrature  sum  gives 

±lij - — - =  -I(r,±/x) 

u  ^  ^  f 

m'=0  r'=c,5  j'  =  i 

+  ^/LjYim'rr' ~ (4.26) 
j  =  l,...,i\r  m  =  0,  ...,M  r^c^s 

4*4.3  Radiative  Transfer  Scattering  Matrix 

It  is  convenient  to  express  the  single  scattering  phase  matrix  P  as  a  Legendre  series 

Ni 

P(cos0)  -  ^X/P/(cos0)  (4.27) 

/=o 

where  Vi  is  the  Vth  order  Legendre  polynomial  and  xi  the  /’th  Legendre  coefficient 
matrix.  For  the  scalar  (unpolarized)  case,  the  Fourier  modes  of  the  phase  function  are 
calculated  from  the  coefficients  of  the  Legendre  series  using  the  addition  theorem  of  associ¬ 
ated  Legendre  functions.  The  rotation  of  the  frame  of  polarization  precludes  that  method 
for  finding  the  Fourier  modes  of  the  scattering  matrix.  The  method  used  here  is  to  perform 
the  polarization  rotation  in  azimuth  space  and  then  Fourier  transform  the  results  to  get 
the  scattering  matrix  for  each  Fourier  azimuth  mode.  This  method  is  similar  to  that  of 
Ishimaru  et  al.  (1984),  except  the  rotation  is  performed  on  the  Stokes  parameters  rather 
than  the  scattering  amplitudes. 

For  each  pair  of  quadrature  angles  fij  and  fiji  (outgoing  and  incoming  directions)  and 
for  a  number  of  azimuth  angles  =  (/>'  —  (/>,  the  scattering  angle  0  is  found  and  the 
Legendre  series  is  summed  for  the  six  unique  elements  of  the  phase  matrix.  Equation  4.17 
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is  used  to  obtain  the  scattering  matrix  M  with  the  desired  polarization  reference  frame. 
The  azimuth  angle  differences  are  at  equally  spaced  angles 

0<7r  Jc 

A<^fc  =  —  =  (4.28) 

where  is  chosen  so  the  highest  frequency  in  A(j)  is  completely  sampled,  or  so 

=  2l^^52(iVi+4)+il+i^  (4  29) 

The  scattering  matrix  M  is  then  Fourier  transformed  to  obtain  the  Fourier  series  repre¬ 
sentation 

M 

-  <l>)  =  -<!))  + smrn{(f>'  -  (f))]  (4.30) 

m—0 

This  matrix  is  a  Fourier  series  in  </>'  —  (/>,  but  it  must  be  transformed  into  a  Fourier  series 
in  (f)^  and  (^,  such  sis 

M  M 

Min -<f>)=  XI  cos m((?!>)  cos m'((/>')  +  M"„,  cos m(<A)  sin m' ((/»') 

m=0  m'=0 

+  smm{<f>)  cosm’ {<!)')  +  smm{(f>)  sinm'((/>')].  (4.31) 

From  angle  addition  formulae  for  cosine  and  sine,  the  coefficients  for  the  double  Fourier 
series  are  related  to  the  single  Fourier  series  coefficients  by 

=  -M-„,  =  (4.32) 

Since  the  scattering  matrix  only  depends  on  the  difference  in  the  azimuth  between  the 
incoming  and  outgoing  angles,  the  azimuth  Fourier  modes  separate,  ie.  M  depends  only 
on  m  rather  than  both  m  and  m! .  This  allows  the  azimuth  modes  to  be  solved  separately, 
reducing  the  computational  burden.  In  the  scalar  radiative  transfer  case  the  scattering 
function  depends  only  on  cosm{(j)^  —  (f>)^  therefore  the  and  terms  are  zero 

and  the  cosine  and  sine  terms  do  not  mix. 

Equation  4.17  has  some  special  symmetries.  The  upper  left  and  lower  right  two  by 
two  blocks  are  even  functions  in  A</),  while  the  upper  right  and  lower  left  blocks  are 
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odd  functions.  Therefore  the  cosine  matrices  have  off-diagonal  blocks  of  zeros  and  the 
sine  matrices  have  diagonal  blocks  of  zeros.  Another  statement  of  this  symmetry  is  that 
negating  (f)  and  0'  results  in  negating  the  off  diagonal  blocks.  This  allows  the  trivial 
calculation  of  the  scattering  matrix  for  tt  <  A(/>  <  27r  from  the  values  for  0  <  A(f>  <  tt. 
Another  symmetry  is  negating  /x  and  /x'  results  in  negating  the  off  diagonal  blocks.  This 
is  due  to  the  invariance  of  the  scattering  matrix  with  rotation  of  the  coordinate  system 
{fi  -/X  and  (j>  — >  M(/x,  /x',  (/>'  -  (^)  =  M(“/x,  -/x', 

For  a  given  azimuth  mode  m  the  cosine  and  sine  scattering  matrices  can  be  combined 


into  a  single  eight-by-eight  scattering  matrix 


'  in 

(4.33) 


The  blocks  of  zeros  in  this  scattering  matrix  allows  a  computationally  efficient  rear¬ 
rangement  of  the  cosine  and  sine  modes  of  Stokes  parameters.  Writing  the  Stokes  radiance 


Mfi 

Mf2 

Mfs 

MU 

MU 

-Mil 

-M|2 
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MU 

-Mil 
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MU 

Mil 

MU 

MU 

MU 

(4.35) 


(4.36) 
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Conventionally  the  solar  azimuth  is  set  to  zero,  so  both  the  solar  and  thermal  sources  of 
radiation  are  even  functions  and  unpolarized,  or  is  the  only  non-zero  element.  Since 
the  F  and  P  radiance  vectors  are  decoupled,  further  scattering  leaves  the  sine  mode 
radiance  vector  zero.  Therefore,  only  the  first  four-by-four  Fourier  scattering  matrix  need 
be  considered.  For  the  azimuthally  symmetric  case  (m  =  0)  the  sine  terms  in  the  scattering 
matrix  are  zero  and  the  U  and  V  Stokes  parameters  are  zero,  in  this  case,  the  two-by-two 
scattering  matrix  may  be  used. 

If  the  solar  azimuth  is  zero  the  direction  of  the  incident  solar  radiation  is  fio  = 
cos  9o,<l>o  =  The  rearranged  Stokes  vector  for  the  solar  ’’pseudo  source”  is  then 


Fo  u  .-r 


(  Mf  1  \ 
Mil 
Mil 


(4.37) 


\  Mil  j 

where  the  scattering  matrix  elements  are  evaluated  as  described  above  for  the  quadrature 
angles  fij  and  the  solar  angle  /Xq. 


The  Matrix  Formation  of  the  Equation  of  Transfer 


Using  the  notation  the  discretized  radiative  transfer  equation  be¬ 


comes. 


dit, 


~  N 

(jJ 


Att 


(4.38) 


j=l 


where  the  plus  and  minus  superscripts  refer  to  the  quadrature  angles  in  the  downward 
(+)  and  upward  (— )  directions  respectively. 

The  notation  can  be  simplified  by  removing  the  subscripts  and  writing  the  equation 
using  matrices  and  vectors.  The  dependencies  on  both  r  and  m  is  taken  to  be  understood. 
The  radiative  transfer  equation  will  be  numerically  solved  for  each  azimuth  mode  m.  The 
Stokes  radiance  vectors  for  each  quadrature  angle  in  a  hemisphere  are  combined  into  one 
long  radiance  vector 


7  = 


^  ii=i  \ 
ii=2 

V  ) 


Ix  = 


(  i]  \ 

9i 

Yi 


(4.39) 
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The  length  of  this  vector  is  4.  x  N.  The  individual  scattering  matrices  for  each  pair  of 
angles  Hj,  fiji  are  combined  to  form  a  scattering  matrix  that  operates  on  a  radiance  vector 
to  produce  another  vector,  which  is  the  the  discrete  equivalent  of  integration  over  all 
angles.  In  matrix  notation  the  radiative  transfer  equation  becomes 

±T>^==-I^  +  C^^I^  +  C^~r  +  d^.  (4.40) 

dr 

The  elements  of  the  matrices  D  and  are 

|D  I 

(4.41) 

where  the  i  and  i'  indices  refer  to  the  Stokes  parameters. 

4.5  Integrating  the  Radiative  Transfer  Equation 
4.5.1  Finite  Diflference 

The  radiative  transfer  equation  (4.40)  is  a  matrix  first  order  ordinary  differential 

equation.  One  method  of  numerically  integrating  the  differential  equations  to  approximate 

the  derivative  (^)  with  a  finite  difference  scheme.  The  simplest  difference  scheme  is 

^  ^  Hr)  -  Hr  +  At)  ^ 

dr  At  '  ^  ’ 

Although  simple,  this  method  is  very  accurate  if  At  is  chosen  small  enough.  Defining 

lo  =  I{t)  and  I\  =  I{t  +  At)  equation  (4.40)  becomes 

/l+  -  /+  =  AtD-1[-/+  +  C++/+  +  C+-/f  +  a+] 

/f  -  /“  =  ArD-^-iT  +  C-+/o^  +  C— Jf  +  a-]  (4.43) 

The  finite  difference  scheme  requires  the  right  hand  side  of  the  equation  be  defined  in 
terms  of  the  discretized  values  lo  and  Ji.  This  is  nearly  the  same  since  Ar  is  chosen  very 
small.  The  last  equation  can  be  arranged  to  express  the  radiation  emerging  from  a  this 
layer  in  term  of  the  incident  radiation. 


=  [1  -  ArD-^i  -  C++)]/+  +  ArD-^C+-If  +  AtD-V+ 
Iq  =  [i  -  AtD-H1-C— )]/f  +  AtD-^C-+/o+  4-  ArD-V- 


(4.44) 

(4.45) 
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AH _  T/q~ 

T i  It 

Figure  4.3:  Illustration  of  the  interaction  principle.  The  and  /f  on  the  left  represent 
incident  radiation.  The  Iq  and  it  on  the  right  represent  the  radiation  leaving  the  medium. 
iJ,T,  and  S  are  the  reflection,  transmission,  and  source  terms  respectively. 

4,5.2  The  Interaction  Principle 

The  interaction  principle  is  a  simple  method  of  expressing  the  linear  interaction  of 
radiation  with  a  medium.  Radiation  emerging  from  any  medium  can  be  expressed  in  term 
of  the  incident  radiation  plus  the  radiation  generated  within  the  medium.  The  matrix 
formulation  for  the  interaction  principle  is 

/+  =  T+I+  +  R+I{  +  S+ 

I-  =  T-I{  +  R-4  +  S-  (4.46) 

where  T  is  the  transmission  matrix,  R  is  the  reflection  matrix,  and  S  is  the  source  vector. 
The  finite  difference  form  of  the  radiative  transfer  equation  (4.40)  is  the  same  form  as  the 
interaction  principle.  The  reflection  and  transmission  matrices  and  the  source  vectors  for 
the  infinitesimal  layer  can  be  related  to  the  local  properties  of  the  medium  by 

r=t  =  {i-ATD-i(i-c±±)] 

R^  =  ArD-^C^^ 

=  ArD-V±.  (4.47) 

Relating  the  radiative  transfer  equation  to  the  interaction  principle  is  normally  termed 
initialization.  There  are  several  initialization  methods,  the  one  used  here  is  sometimes 
called  the  infinitesimal  generator  initialization. 
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The  Adding  Algorithm 

The  previous  section  have  shown  how  the  coefficients  of  the  radiative  transfer  equation 
relate  to  the  reflection  matrices,  transmission  matrices,  and  source  vectors  for  infinitesimal 
layers.  Now  an  integration  procedure  is  needed  to  relate  the  properties  of  the  infinitesimal 
layers  to  the  properties  of  an  atmosphere  composed  of  many  infinitesimal  layers.  With 
the  reflection  and  transmission  matrices  and  the  source  vectors  for  the  entire  medium,  the 
interaction  principle  can  be  applied  to  the  incident  radiation  at  the  boundaries  to  calculate 
the  outgoing  radiation.  Radiation  fields  inside  the  medium  may  also  be  calculated  from 
the  incident  radiation  and  the  matrices  and  source  vectors  for  the  medium  above  and 
below  the  level  of  interest. 


Figure  4.4:  Illustration  of  adding  principle.  The  properties  of  layers  1  and  2  are  combined 
to  make  the  reflection,  transmission,  and  source  term  for  the  new  layer. 

The  interaction  principle  is  now  used  to  combine  (add)  two  layers.  The  reflection 
and  transmission  matrices  and  the  source  vector  for  this  new  double  thick  layer  can  be 
expressed  in  terms  of  the  matrices  and  vectors  of  the  individual  layers.  Starting  with 
the  interaction  principle  for  the  two  layers,  Figure  4.4,  the  radiances  at  the  interface  are 
eliminated  and  the  equations  rearranged  to  form  the  interaction  principle  for  the  new 
layer,  resulting  in 

R^  =  Rt  +  T^r+RfT^,  +  T^V-R^T^ 

=  T^T+T^,  Tf  =  Tf  r-r2" 

5^  =  5^  +  r2+r+(5J^  +  Rts^),  Sj  =  5f  +  Tf  r-(52“  +  R^St) 

r+  =  [1  -  r-  =  [1  -  R^Rtr\  (4.48) 
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By  similar  manipulation,  the  internal  radiances  {I^-,  Ii  )  may  be  expressed  in  term  of 
the  radiance  incident  upon  the  two  layers  as 

1+  =  +  T+I+  +  R+S^  +  5+1 

/f  =  r- [i22  Ti+ Jo+  +  T^l2  +  5i+  +  52  ].  (4.49) 

Doubling  Algorithm 

The  radiative  properties  for  a  finite  layer  could  be  calculated  by  building  on  one  layer 
at  a  time  using  the  adding  algorithm,  however,  the  doubling  algorithm  provides  a  much 
faster  method  of  building  up  a  thick  layer  from  many  identical  thin  layers.  If  all  the  thin 
layers  have  an  optical  depth  of  At,  then  combining  the  first  two  layers  yields  a  layer  with 
optical  depth  2Ar.  Since  all  the  thin  layers  are  identical,  all  the  2Ar  layers  will  be  the 
same.  Two  2Ar  layers  are  now  added  together  to  form  a  4At  thick  layer.  After  N  of 
these  doubling  steps  the  layer  is  built  up  to  an  optical  depth  of  2^  Ar.  Using  this  method 
more  than  one  million  thin  layers  can  be  combined  in  20  steps.  The  doubling  formulae  for 
the  reflection  and  transmission  matrices  are 

^2N  ”  ^n'^N  ’  ^2N  “ 

T2V  =  T+r+r  +  N+,  T2I,  = 

T+  =  [l-R+R],]-\  T-  =  [l-R-^R+]-\  (4.50) 

The  doubling,  adding,  and  finite  generator  initialization  algorithms  used  here  were 
developed  by  Grant  and  Hunt  (1969).  The  doubling  method  so  far  requires  all  the  layers 
be  uniform. 

The  solar  pseudo-source  has  an  exponential  dependence  with  optical  depth.  Wis- 
combe  (1976)  extended  the  doubling  method  to  incorporate  these  sources.  With  exponen¬ 
tial  sources  the  source  at  a  deeper  optical  depth  is  the  same  as  the  source  at  a  shallower 
optical  depth  except  for  a  multiplicative  factor.  This  introduces  a  factor  into  the  formula 
for  adding  the  sources  of  two  layers 


^2N  ”  'y^ ^iv'^ArT^) 
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52lv  =  +  R-^S],) 

7  =  exp{-AT/fjio).  (4.51) 

4.5.3  The  Boundaries 

The  effect  of  the  surface  may  also  be  calculated  using  the  interaction  principle.  The 
ground  is  treated  as  a  layer  with  transmission  of  unity,  no  source,  and  a  reflection.  The 
ground  is  assumed  to  be  a  Lambertian  surface,  and  therefore  emits  and  reflects  equally  in 
all  directions.  The  reflected  radiance  is  the  incident  flux  times  the  ground  albedo  divided 
by  It.  Reflected  radiance  is  assumed  to  be  unpolarized.  Radiation  emitted  from  the  ground 
is  the  incident  radiance  on  the  lower  boundary.  The  internal  radiance  algorithm  is  used 
to  calculate  the  downwelling  radiance  below  the  atmosphere,  but  above  the  ground 

/-  =  r+/+' 

if  =  T+I+  =  </-  +  S+ 

r+  =  [1  -  (4.52) 

where  a  refers  to  the  total  atmosphere  layer  and  g  refers  to  the  ground.  The  downwelling 
radiation  from  the  atmosphere  is  what  the  radiation  would  be  with  no  ground  multi¬ 
plied  by  the  multiple  reflection  factor  . 

The  radiance  upwelling  from  the  top  of  the  atmosphere  is  found  by  adding  the  atmo¬ 
sphere  layer  to  the  ground  layer 

+  Trp  Ig  =  Srq 

Rt  =  '^a^^Rg'^a 

TT=T-r- 


Srp  ^  H-  Tq  r  RgS^ 
r-  =  [1  -  RgRt]-\ 


(4.53) 
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4.6  Testing  the  Polarized  Radiative  Transfer  Model 

A  convenient  test  of  the  model  is  to  compare  the  results  to  the  tables  by  Coulson  et  al. 
(1960).  The  Coulson  tables  list  the  outgoing  radiation  for  a  conservative  sunlit  Rayleigh 
atmosphere.  Both  the  upwelling  and  downwelling  Stokes  radiance  vectors  (J,  Q,i7)  are 
provided  in  terms  of  different  optical  depths,  solar  angles,  and  surface  albedos.  Compar¬ 
isons  were  made  with  two  cases  of  varying  optical  depth  and  solar  angles.  The  radiative 


Table  4.1:  Comparison  of  model  results  with  Coulson  tables  for  a  homogeneous  Rayleigh 
atmosphere.  The  upwelling  radiance  at  azimuth  angle  of  90°  for  optical  depth  of  1,  solar 
zenith  angle  of  0.8,  and  surface  albedo  of  0.25.  Coulson  defines  Q  with  a  sign  opposite  of 
that  used  in  this  model. 


Coulson 

Model 

I 

Q 

U 

I 

Q 

U 

0.06 

0.39887 

0.05099 

0.24758 

0.395222 

-0.05162 

0.246095 

0.16 

0.40894 

0.03988 

0.23375 

0.406471 

-0.04037 

0.232832 

0.28 

0.40482 

0.02766 

0.20918 

0.403094 

-0.02801 

0.208589 

0.4 

0.3938 

0.0157 

0.17114 

0.392471 

-0.01596 

0.180708 

0.64 

0.37248 

-0.00774 

0.12476 

0.371605 

0.007621 

0.124514 

0.84 

0.36147 

-0.02681 

0.0759 

0.360826 

0.026772 

0.075762 

0.96 

0.35776 

-0.03808 

0.03609 

0.357222 

0.038087 

0.036026 

1 

0.35694 

-0.04181 

0 

0.356435 

0.041833 

0 

transfer  model  was  run  with  16  angles  using  the  extra  angle  quadrature  mode  (8  angles 
machine  generated  and  8  user  supplied)  and  azimuth  modes  up  to  m  =  2.  The  upwelling 
and  downwelling  radiances  were  compared  at  azimuth  angles  of  0,  90,  and  180°  (U  is  zero 
at  0  and  180°).  Table  4.1  compares  the  upwelling  radiance  at  azimuth  angle  of  90°  for  a 
homogeneous  Rayleigh  atmosphere  with  optical  depth  of  1,  solar  zenith  angle  of  0.8,  and 
surface  albedo  of  0.25.  Note:  Coulson  define  Q  with  a  sign  opposite  of  that  used  in  this 
model.  Table  4.2  compares  the  downwelling  radiance  at  azimuth  angle  of  180°  for  optical 
depth  of  1,  solar  zenith  angle  of  0.6,  and  surface  albedo  of  0.  Table  4.3  summarizes  the 
two  cases,  the  average  and  maximum  differences  between  the  model  and  Coulson  tables 
are  shown. 
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Table  4.2:  Comparison  of  model  results  with  Coulson  tables  for  a  homogeneous  Rayleigh 
atmosphere.  The  upwelling  radiance  at  azimuth  angle  of  180®  for  optical  depth  of  .15, 
solar  zenith  angle  of  0.6,  and  surface  albedo  of  0.  U  is  zero  at  0  and  180®. 


■1 

Coulson 

Model 

I 

Q 

I 

Q 

0.06 

0.29264 

0.05675 

0.288185 

-0.05806 

0.16 

0.18709 

0.05065 

0.184099 

-0.05142 

0.28 

0.11837 

0.04489 

0.116442 

-0.04528 

0.4 

0.08206 

0.0414 

0.080709 

-0.04158 

0.64 

0.04656 

0.03616 

0.04583 

-0.03608 

0.84 

0.03438 

0.03041 

0.033906 

-0.03027 

0.96 

0.03347 

0.02387 

0.033064 

-0.02373 

1 

0.03859 

0.01662 

0.038152 

-0.01654 

Table  4.3:  Summary  of  differences  between  the  model  and  Coulson  tables.  The  average 
and  maximum  absolute  differences  of  the  radiances  over  the  eight  zenith  angles  at  azimuths 
of  0,  90,  180®  and  both  upwelling  and  downwelling  directions. 


Optical 

Solar 

Surface 

Average  Error 

Maximum  Error 

Depth 

Albedo 

I 

Q 

U 

I 

Q 

U 

1 

0.8 

0.25 

0.00108 

0.00014 

0.00024 

0.00399 

0.00063 

0.00149 

0.15 

0.6 

0 

0.00174 

0.00057 

0.00021 

0.00519 

0.00211 

0.00193 
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4.7  Model  Results 
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Figure  4.5:  Scattering  pattern  for  a  Rayleigh  particle.  Incident  radiation  is  unpolarized. 
//  is  the  component  parallel  to  the  scattering  plane  and  //  is  the  component  perpendicular 
to  the  scattering  plane.  The  pattern  is  rotationally  symmetric  about  the  direction  of  the 
incident  radiation. 


For  unpolarized  incident  radiation,  the  angular  distribution  of  Rayleigh  scattering  is 
shown  in  Figure  4.5.  The  distribution  has  rotational  symmetry  around  both  0  =  0—180° 
and  0  =  90°.  Radiation  scattered  directly  forward  or  backward  remains  totally  unpolar¬ 
ized.  Radiation  scattered  at  0  =  90°  has  half  the  intensity  of  forward  and  back  scattered 
radiation,  but  is  totally  polarized.  The  definition  of  the  Rubenson  degree  of  polarization 
(Coulson,  1988)  comes  from  this  relationship 


P=: 


Ir-h 

Ir  +  II 


(4.54) 


This  definition  of  linear  polarization  {V  =  0)  is  convenient  when  in  the  plane  of  the  sun’s 
vertical,  since  it  is  eetsy  to  calculate  and  allows  for  both  positive  and  negative  polariza¬ 
tion.  However,  at  directions  that  are  not  in  the  plane  of  the  sun’s  vertical  the  Rubenson 
definition  gives  erroneous  results  because  Ir  and  //  with  respect  to  the  viewing  direction 
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are  not  the  same  as  Ir  and  //  with  respect  to  the  scattering  plane.  At  directions  not  in 
the  sun’s  vertical  the  state  of  linear  polarization  in  given  by 


Ir  +  h 


(4.55) 


with 


U  —  {Ir  —  Ii)  tan  2x 


(4.56) 


where  x  the  angle  the  plane  of  polarization  taken  with  respect  to  the  vertical  direction. 
The  degree  of  polarization  is  normally  positive  due  to  the  square  root  relation.  Only  with 
i7  =  0  does  a  negative  sign  of  P  have  meaning.  Since  i7  =  0  in  the  suns  vertical,  P  can 
become  negative  the  Rubenson  definition  applies.  In  the  most  general  sense  polarization  is 
not  linear  {V  7^  0),  using  the  (I,  Q,  V)  representation  the  degree  of  polarization  becomes 

+  (4.57) 

This  representation  of  polarization  does  not  allow  negative  polarization,  but  allows  calcu¬ 
lation  of  the  degree  of  polarization  at  all  angles.  It  is  the  relationship  used  to  calculate 
polarization  for  the  sky  hemisphere. 

Figure  4.6  shows  the  radiation  intensity  for  a  pure  Rayleigh  atmosphere.  The  solar 
zenith  angle  is  78.5°,  optical  depth  is  .05  and  surface  albedo  is  zero.  This  is  the  familiar 
Rayleigh  pattern  for  atmospheric  radiation  intensity  with  it’s  nearly  concentric  contours 
and  the  brightening  near  the  horizon.  Figure  4.7  shows  the  sky  polarization  for  the  same 
solar  angle,  optical  depth,  and  surface  albedo.  Due  to  the  low  sun  angle  and  relatively 
small  optical  depth  there  is  little  multiple  scatter  and  a  95%  polarization  band  occurs 
close  to  the  90°  azimuth. 

Figures  4.8  and  4.9  show  the  sky  polarizations  for  solar  zenith  angles  of  66.4°  and 
23.1°  respectively.  For  these  two  causes,  the  optical  depth  wa5  increased  to  .15,  but  surface 
albedo  is  held  at  zero.  In  both  cases  we  see  the  effect  of  scattering  angle  the  maximum 
polarization  band  remains  close  to  90°  from  the  sun.  The  maximum  polarization  band 
shifts  across  the  sky  and  the  contours  become  more  concentric  with  decreasing  solar  zenith 
angle.  If  the  sun  were  at  the  zenith,  the  maximum  polarization  would  be  at  the  horizon 
and  the  polarization  contours  would  be  symmetric  about  the  solar  zenith. 
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Figure  4.6:  Distribution  of  the  transmitted  intensity  over  the  hemisphere  of  the  sky.  The 
solar  position  is  azimuth  =  0  and  zenith  angle  =  78.5°,  optical  depth  is  .05  and  surface 
albedo  is  zero.  Only  half  the  hemisphere  is  shown,  the  other  half  being  the  mirror  image 
of  this  one.  The  plane  of  the  sun’s  vertical  is  represented  by  the  straight  bottom  line,  with 
0  being  the  zenith.  The  outside  semicircle  represents  the  horizon. 

Figure  4.10  shows  the  hemispheric  radiation  intensity  with  the  introduction  of  haze 
into  the  atmosphere.  Here  the  solar  zenith  angle  is  again  23.1°  and  the  surface  albedo 
is  zero.  The  optical  depth  is  also  .15,  however  .10  is  due  to  Rayleigh  scattering  and  .05 
is  due  to  the  haze.  The  intensity  close  to  the  sun  sharply  increases  due  to  the  strong 
forward  scatter  for  the  haze.  Intensity  at  all  other  angles  also  increases  due  to  increased 
multiple  scattering.  Figure  4.11  shows  the  corresponding  graph  for  the  sky  polarization. 
The  band  of  maximum  polarization  does  not  change  location  much,  although  there  is  a 
decrease  in  the  maximum  intensity.  There  is  also  a  larger  area  of  low  polarization  near 
the  solar  position. 

To  closely  examine  the  effects  of  the  parameters  that  affect  the  sky  radiation  intensity 
and  polarization,  the  directions  ^  =  0  =  180  were  used.  Figure  4.12  shows 

the  effect  of  increasing  Rayleigh  optical  depth  on  the  sky  light  intensity  and  polarization. 
For  this  exercise,  the  solar  zenith  angle  was  held  constant  at  50°  (solar  elevation  =  40°) 
and  the  surface  albedo  was  zero.  The  sky  intensity  overall  shows  a  minimum  at  high 
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Figure  4.7:  Distribution  of  the  degree  of  polarization  over  the  hemisphere  of  the  sky.  The 
solar  position  is  azimuth  =  0  and  zenith  angle  =  78.5°,  optical  depth  is  .05  and  surface 
albedo  is  zero. 

elevation  angles  and  a  maximum  at  low  elevation  angles  {9  <  40°  and  9  >  140°).  This 
is  a  result  of  a  plane  parallel  atmosphere.  At  high  elevation  angles,  the  path  length 
through  the  atmosphere  is  short  and  not  much  radiation  is  scattered  from  the  direct  solar 
beam  into  the  viewing  direction.  At  low  elevation  angles,  the  path  length  through  the 
atmosphere  increases  and  more  radiation  is  scattered  from  the  direct  solar  beam  into 
the  viewing  direction.  As  optical  depth  increases,  sky  radiation  intensity  increases  at  all 
elevations,  since  more  radiation  is  scattered  out  of  the  direct  solar  beam  with  increased 
optical  depth.  The  increase  is  the  most  significant  in  the  sun’s  hemisphere  and  close  to 
the  horizon  since  Rayleigh  scattering  is  stronger  in  the  forward  and  backward  directions. 
The  effects  on  polarization  are  more  complicated.  For  a  purely  single  scattering  Rayleigh 
case,  there  would  be  totally  unpolarized  radiation  (zero  polarization)  at  the  solar  elevation 
and  totally  polarized  radiation  (polarization  =  1)  at  90°  above  the  solar  angle.  However, 
the  atmosphere  is  not  single  scattering.  Multiple  scattering  has  the  effect  of  converting 
perpendicular  polarized  radiation  into  parallel  polarized  radiation.  Therefore  at  elevations 
close  to  the  solar  position,  multiple  scattering  will  cause  the  observed  polarization  to 
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Figure  4.8:  Distribution  of  the  degree  of  polarization  over  the  hemisphere  of  the  sky.  The 
solar  position  is  azimuth  =  0  and  zenith  angle  =  66.4°,  optical  depth  is  .15  and  surface 
albedo  is  zero. 

become  negative.  At  elevations  close  to  90°  from  the  solar  position,  multiple  scattering  will 
reduce  the  observed  polarization  below  the  maximum  of  1  for  single  scatter.  Figure  4.12 
shows  that  as  the  optical  depth  increases,  and  therefore  multiple  scatter  increases,  the 
polarization  near  the  sun  becomes  more  negative,  and  the  polarization  decreases  at  90° 
from  the  sun. 

Figure  4.13  shows  the  effects  of  surface  reflection.  As  the  surface  albedo  increases, 
there  is  an  almost  uniform  increase  in  intensity  for  all  elevation  angles,  this  increase  is 
of  less  magnitude  than  that  for  increasing  the  optical  depth.  The  intensity  increase  is 
a  result  of  radiation  being  reflected  from  the  ground  and  then  being  scattered  into  the 
viewing  direction.  The  polarization  is  effected  the  most  at  angles  90°  from  the  solar 
position.  Increasing  surface  reflection  decreases  the  polarization  maximum  significantly. 
At  elevations  close  to  the  solar  position  there  is  very  little  change  in  polarization.  This  is 
due  to  the  surface  acting  bs  a  lambertian  reflector.  Radiation  reflected  from  the  surface 
back  into  the  atmosphere  will  be  unpolarized,  and  induce  a  depolarization  effect  at  all 
angles.  There  will  be  a  large  effect  90°  from  the  sun  at  the  polarization  max,  but  since 
the  polarization  near  the  sun  is  already  close  to  zero  the  effect  here  will  be  less  noticeable. 
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Figure  4.9:  Distribution  of  the  degree  of  polarization  over  the  hemisphere  of  the  sky.  The 
solar  position  is  azimuth  =  0  and  zenith  angle  =  23.1"^,  optical  depth  is  .15  and  surface 
albedo  is  zero. 

Figure  4.14  shows  the  effect  solar  zenith  angle  has  on  the  sky  light  intensity  and 
polarization.  As  the  solar  zenith  angle  decreases  (solar  position  closer  to  90°)  the  intensity 
pattern  becomes  more  symmetric.  As  the  sun  changes  position,  the  areas  of  maximum  and 
minimum  polarization  track  the  solar  position.  As  the  solar  position  moves  higher  into  the 
sky  the  area  of  negative  polarization  near  the  solar  position  becomes  less  negative,  due  to 
optical  path  decreasing,  and  therefore  less  multiple  scattering.  The  opposite  happens  at 
90°  from  the  solar  position.  As  the  solar  position  rises  into  the  sky  the  area  of  maximum 
polarization  now  has  an  increasing  optical  path  and  the  resulting  increase  in  multiple 
scattering  causes  the  observed  maximum  polarization  to  decrease  in  magnitude. 

Figure  4.15  shows  the  effect  of  adding  haze  to  the  bottom  kilometer  in  the  atmosphere. 
For  this  test  the  Rayleigh  scatter  Wcts  held  at  .10  and  the  optical  depth  due  to  the  haze  was 
increased  from  .01  up  to  .10.  The  solar  zenith  angle  was  50°  and  the  surface  albedo  was 
zero.  A  Deirmendjian  Haze-L  distribution  (Deirmendjian,  1969)  was  used  with  a  complex 
index  of  refraction  m  =  1.50  —  .OOOi  at  A  =  520nm.  The  scattering  phase  function  was 
determined  by  mie  calculations,  particles  are  assumed  to  be  randomly  oriented.  The 
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Figure  4.10:  Distribution  of  radiation  intensity  over  the  hemisphere  of  the  sky.  The  solar 
position  is  azimuth  =  0  and  zenith  angle  =  23.1°,  Rayleigh  optical  depth  is  .10,  the  haze 
optical  depth  is  .05,  and  surface  albedo  is  zero. 

intensity  phase  function  is  shown  in  Figure  4.16  with  the  phase  function  for  Rayleigh 
scattering.  We  see  the  haze  has  very  strong  forward  scattering  peak  from  0°  through  40° 
and  the  backscatter  is  less  than  that  for  Rayleigh  scattering.  The  strong  forward  scattering 
of  the  haze  shows  up  readily  in  Figure  4.15  as  the  narrow  peak  in  intensity  near  the  sun. 
Away  from  the  solar  position,  the  increa.se  in  intensity  falls  off  rapidly.  Haze  de-polarizes 
the  the  Rayleigh  atmosphere  at  all  angles.  In  the  vicinity  of  the  solar  position,  multiple 
Rayleigh  scattering  normally  introduces  a  negative  polarization,  but  haze  over  shadows 
this  effect  and  the  resultant  polarization  is  close  to  zero.  At  90°  from  the  solar  position 
there  is  a  larger  de-polarization  as  the  optical  depth  of  the  haze  increases.  The  point  of 
maximum  polarization  also  moves  further  from  the  solar  position  as  the  haze  optical  depth 
increases.  This  is  a  result  of  the  haze  forward  scattering  depolarizing  the  area  closer  to 
the  sun  to  a  greater  degree  than  those  area  farther  away  from  the  sun. 


ZENITH  ANGLE  (deg) 

Figure  4.11:  Distribution  of  the  degree  of  polarization  over  the  hemisphere  of  the  sky.  The 
solar  position  is  azimuth  =  0  and  zenith  angle  =  23.1°,  Rayleigh  optical  depth  is  .10,  the 
haze  optical  depth  is  .05,  and  surface  albedo  is  zero. 
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Figure  4.13;  Effect  of  changing  surface  albedo  on  sky  radiation  intensity  and  polarization. 
Solar  position  is  held  at  40°  and  the  optical  depth  is  .15. 
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Figure  4.15:  Effects  of  introducing  haze  into  the  atmosphere.  The  Rayleigh  portion  of 
optical  depth  is  held  constant  at  .10  while  the  optical  depth  of  haze  increased  from  .01  to 
.10.  Solar  zenith  angle  is  50°  and  surface  albedo  is  zero. 


Chapter  5 


SKY  POLARIZATION  EXPERIMENT 

In  this  experiment,  the  sky  radiation  intensity  is  measured  with  the  SSP.  Using  the 
calibration  coefficients  described  earlier,  the  intensity  values  are  used  to  calculate  sky 
polarization.  The  resultant  polarization  patterns  are  compared  to  the  modeled  values. 

5.1  Measuring  the  Solar  Radiation 

The  SSP  was  mounted  in  the  CSU  Smart  Solar  Tracker  as  shown  in  Figure  5.1.  The 


Figure  5.1:  CSU  Smart  Solar  Tracker. 
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solar  tracker  is  a  servo-driven  two-axis  mount  that  is  controlled  by  a  PC  based  motor 
controller  card.  The  position  of  the  servo  motors  is  determined  using  optical  encoders, 
which  have  a  resolution  of  600  steps/degree  of  rotation.  Solar  position  is  calculated  using 
the  positioning  algorithm  from  ’’Astronomical  Algorithms”  (Meeus,  1991).  Two  modes  of 
the  tracking  software  were  used.  In  the  first  mode  the  tracker  centers  itself  on  the  sun  and 
continually  adjusts  its  position  to  follow  the  sun’s  movement.  The  second  mode  allows 
the  user  to  stop  the  tracking  routine  and  have  the  tracker  move  a  given  number  of  degrees 
from  the  solar  position. 

The  procedure  for  this  exercise  included: 

1.  Set  the  SSP  in  the  solar  tracker.  Align  the  SSP  center  line  in  the  vertical.  Cover 
the  SSP  lenses  and  center  the  solar  tracker  on  the  sun. 

2.  Once  the  solar  tracker  has  located  the  sun,  stop  the  tracking  software.  Use  the 
controller  to  manually  move  the  solar  tracker  a  given  number  of  degrees  from  the  solar 
position.  These  moves  are  always  made  in  the  elevation  of  the  tracker.  Uncover  the  SSP 
lenses  and  record  the  sky  radiation  intensity. 

3.  Cover  the  lenses  and  have  the  tracking  software  re-acquire  the  sun. 

4.  Repeat  Steps  2  and  3. 

Intensity  measurements  were  always  made  in  the  vertical  plane  of  the  sun.  This 
simplifies  several  tasks:  the  solar  tracker  only  has  to  be  adjusted  in  elevation,  calculation 
of  the  scattering  angle  is  trivial  (it  is  simply  the  difference  between  the  solar  elevation 
and  the  elevation  the  solar  tracker  is  pointed),  polarization  can  be  determined  using  the 
Rubenson  formula,  and  both  the  maximum  and  minimum  sky  polarizations  will  be  sampled 
during  each  scan. 

Measurements  were  made  at  5^  intervals  in  elevation  both  above  and  below  the  solar 
position.  At  each  position  data  were  recorded  for  approximately  60  seconds.  This  provided 
more  than  100  scans  of  data  to  be  averaged  together.  Additional  time  was  requii^ed  to 
re-acquire  the  solar  position  before  each  reading.  To  scan  the  entire  hemisphere  required 
about  2  hours.  Three  scans  were  performed  on  24  February,  the  first  and  third  scans 
centered  around  2^  hours  before  and  after  solar  noon  with  an  average  solar  elevation  of 
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30^  and  the  second  scan  wcis  centered  around  solar  noon  with  the  solar  elevation  of  40°. 
Data  from  the  second  scan  will  be  used  for  illustration  purposes. 

During  the  same  time  the  sky  radiation  was  being  measured,  measurements  of  optical 
thickness  were  made  with  the  University  of  Arizona,  Custom  Filter- Wheel  Solar  Radiome¬ 
ter.  The  radiometer  uses  ten  narrow  band  filters  from  370  nm  through  1030  nm.  Data 
from  the  radiometer  are  plotted  on  a  langley  diagram  to  determine  the  optical  thickness. 
Table  5.1  lists  the  measured  optical  thickness  for  each  of  the  radiometer  channels. 


Table  5.1:  The  University  of  Arizona,  Custom  Filter-Wheel  Solar  Radiometer  central 
wavelengths  and  the  measured  optical  depths. 


Channel 

Wavelength 
Center  (nm) 

Optical 

Thickness 

1 

.431 

2 

.317 

3 

441 

.221 

4 

.133 

5 

.097 

6 

.057 

7 

779 

.027 

8 

.018 

9 

939 

.193 

10 

.012 

5,2  Solar  Data 

Figure  5.2  illustrates  an  interesting  point  about  the  calibration  lamp.  The  plot  shows 
the  raw  data  signal  received  from  the  calibration  source  at  a  gain  setting  of  25.75.  and 
the  raw  power  received  from  the  sky  when  the  SSP  is  10°  from  the  solar  position  at 
the  same  gain  setting.  Note  how  the  received  power  from  the  calibration  lamp  is  much 
lower  than  the  received  power  from  the  sky  light  at  wavelengths  shorter  than  700  nm. 
Thus  for  the  shorter  wavelengths,  error  may  be  introduced  into  the  calibration  coefficients 
due  to  instrument  noise,  which  is  the  spectral  region  where  the  atmospheric  signal  is  the 
strongest.  The  only  way  to  overcome  this  is  to  find  a  different  type  of  calibration  source 
that  emits  stronger  at  the  shorter  wavelengths. 
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Figure  5.2:  Raw  data  signals  received  by  the  SSP  for  the  calibration  lamp  and  from  the  the 
sky  with  the  SSP  10°  from  the  solar  position.  The  gain  setting  for  both  meaisurements  was 
25.75.  Note  the  low  signal  intensity  from  the  calibration  lamp  at  the  shorter  wavelengths. 

Before  we  examine  the  data  results  from  the  solar  scan,  it  is  interesting  to  once  again 
look  at  the  raw  data.  Figure  3.15,  shown  earlier,  Wcis  for  a  solar  elevation  only  10°  from 
the  solar  position.  Figure  5.3  illustrates  the  same  data  for  the  solar  elevation  90°  from  the 
solar  position.  For  this  scan,  the  gain  setting  on  channels  1  and  2  is  25.75  and  the  gain 
setting  for  channels  3  and  4  is  34.00.  The  higher  gain  setting  for  channels  2-4  introduces  a 
significant  amplification  of  the  DC  voltage  offset,  moving  the  data  points  down  the  graph. 
More  importantly,  the  signal  from  the  sky  is  so  weak  the  spectral  signal  is  difficult  to 
see  visually.  Also,  the  broadband  signal  no  longer  saturates  the  detectors.  This  weak 
signal  is  the  main  reason  60  seconds  (more  than  100  scans)  of  data  were  gathered  for  each 
measurement. 

Figure  5.4  shows  the  measured  spectral  radiance  for  an  entire  scan  through  the  sun’s 
vertical  plane.  This  graph  is  plotted  backward  from  the  previous  graphs,  with  longer 
wavelengths  at  the  front.  For  this  scan,  the  solar  elevation  was  40°  (solar  zenith  angle 
of  50°).  There  is  a  increase  in  radiance  at  the  shorter  wavelengths,  as  is  expected  when 
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Figure  5.3:  SSP  raw  data  at  90"^  from  solar  position.  Gain  setting  on  channels  1  and  2  is 
25.75.  Gain  setting  for  channels  3  and  4  is  34.00.  The  zero  signal  point  is  at  the  top  of 
the  graph. 

measuring  the  solar  spectrum.  The  SSP  also  detects  the  increase  in  radiance  as  it  views 
closer  to  the  horizon,  especially  near  165^  elevation.  Also  note  that  at  5°  elevation,  the 
radiance  falls  off  sharply,  since  the  SSP  was  viewing  a  building.  The  elevations  near  the 
solar  position  are  also  interesting.  Radiance  increases  at  angles  close  to  the  sun.  This  is 
due  to  the  forward  scattering  of  haze.  The  haze  contribution  is  in  a  very  narrow  band 
band  around  the  solar  position.  There  is  even  a  trough  of  low  radiance  values  at  15  - 
20°  elevation  where  the  forward  scattering  of  the  haze  has  not  overshadowed  the  normal 
Rayleigh  radiation  decrease  for  this  area. 
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Figure  5.4:  Spectral  intensity  versus  wavelength  and  elevation  angle  for  24  February.  The 
solar  elevation  was  40°. 


Figure  5.5  shows  the  spectral  intensity  with  its  associated  error  at  60°  and  130° 
elevation.  Several  interesting  items  are  evident.  The  spectral  intensity  is  much  weaker  at 
130°.  This  effect  is  stronger  as  wavelength  increases,  since  Rayleigh  optical  depth  decreases 
with  increasing  wavelength.  There  is  very  little  radiance  in  the  NIR  wavelengths  at  130°. 
The  radiance  values  calculated  for  CVF  sections  1  and  2  in  the  region  of  overlap,  between 
660  nm  and  740  nm,  give  almost  identical  results.  This  is  an  indication  of  the  accuracy 
of  the  relative  calibration.  The  oxygen  absorption  band  at  760  nm  is  detected.  The  water 
vapor  absorption  band  at  940  nm  is  detected  by  the  958  nm  reading.  This  data  shown 
here  is  from  channel  3  for  which  958  nm  is  the  closest  central  wavelength,  however,  the 
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Figure  5.5:  Spectral  intensity  at  60°  and  130°  elevation.  The  solar  elevation  was  40°.  Note 
the  radiance  values  calculated  from  the  different  filter  wheels  between  660  nm  and  740  nm 
give  almost  identical  results.  Also,  the  oxygen  absorption  band  at  760  nm  and  the  water 
vapor  absorption  band  at  940  nm  are  detected. 

HBW  is  wide  enough  that  information  from  the  water  vapor  absorption  band  is  measured. 
As  discussed  earlier,  radiance  errors  are  largest  for  the  shorter  wavelength  due  to  the  low 
emittance  of  the  calibration  lamp  in  this  spectral  region.  In  general,  CVF  section  2  has 
smaller  errors  since  its  transmission  is  higher.  The  errors  for  the  NIR  wavelengths  of 
the  130°  elevation  scan  are  significantly  larger  than  for  the  same  wavelengths  in  the  60° 
elevation  scan  due  to  the  low  signal  strength  at  130°,  thus  allowing  instrument  noise  to 
be  a  larger  fraction  of  the  measurement. 

Figure  5.6  shows  the  spectral  polarization  for  24  February  versus  the  elevation  angle 
and  wavelength.  The  polarization  data  are  noisy,  but  the  maximum  polarization  for 
all  wavelengths  stands  out  at  approximately  130°  elevation,  which  is  90°  from  the  solar 
position.  There  is  also  a  band  of  minimum,  even  slightly  negative,  polarization  near 


0 

Figure  5.6:  Spectral  polarization  versus  wavelength  and  elevation  angle.  The  solar  eleva¬ 
tion  was  40°. 

the  solar  position.  The  polarization  at  all  angles  shows  very  little  dependance  on  the 
wavelength  and  therefore  the  optical  depth. 

Since  all  measurements  are  made  in  the  plane  of  the  sun’s  vertical,  the  polarization 
can  be  calculated  with  Rubenson’s  formula,  or 

p  Ir-h 

Ir  +  h' 

The  error  in  calculating  the  polarization  depends  on  two  quantities,  both  Ir  and  //.  As¬ 
suming  the  errors  are  independent  of  each  other  the  final  error  is  calculated  by 


r 
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which  reduces  to 

(5-2) 

Figure  5.7  shows  the  spectral  polarization  and  its  associated  error  at  60°  and  130° 
elevation.  The  errors  are  substantial  at  all  wavelengths.  At  the  shorter  wavelength  instru¬ 
ment  noise  introduces  errors  into  the  calibration  coefficients.  At  the  longer  wavelengths, 
the  low  sky  light  intensity  due  to  less  scattering  allows  instrument  noise  to  again  intro¬ 
duce  error.  The  net  result  is  the  error  for  the  polarization  measurements  is  large  for  all 
wavelengths.  This  again  is  justification  for  the  long  dwell  time  for  the  sky  intensity  mea¬ 
surements.  Also,  since  intensity  measurements  from  both  channels  2  and  4  are  required  to 
determine  the  polarization,  only  data  points  where  the  central  wavelengths  of  these  two 
channels  are  within  1%  of  each  other  is  the  polarization  calculated.  This  eliminates  some 
of  the  NIR  data  points  shown  in  previous  plots. 

To  evaluate  the  error,  the  distribution  of  individual  radiance  values  at  130°  elevation 
is  shown  in  Figure  5.8  for  all  three  NFOV  channels.  The  range  between  the  maximum 
and  minimum  radiance  value  for  each  channel  was  divided  into  20  equal  bins  and  the 
number  of  readings  within  each  bin  is  plotted.  The  radiance  distribution  for  each  channel 
is  close  to  a  gaussian  distribution.  Skewness  for  a  gaussian  curve  is  zero.  The  skewness 
for  channels  2-4  was  .19,  .14,  and  .1  respectively.  This  indicates  the  majority  of  the 
error  is  associated  with  instrument  noise  and  the  distribution  average  is  a  representative 
value.  The  data  at  130°  elevation  has  a  signal  to  noise  ration  of  4:1  for  the  radiance  and 
perpendicular  polarization  channels  and  as  low  as  1:1  for  the  parallel  polarization  channel. 

5,3  Model  Comparisons 

To  run  the  model,  several  parameters  must  be  determined.  These  include  the  solar 
zenith  angle,  Rayleigh  and  haze  optical  depths,  surface  albedo,  and  solar  spectral  flux  at 
the  top  of  the  atmosphere. 

The  solar  zenith  angle  for  the  second  scan  was  50°.  During  the  entire  scan  this  only 
changed  by  less  than  For  scans  one  and  three,  the  solar  zenith  angle  varied  from 
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56°  through  68°.  For  these  scans,  an  average  solar  zenith  angle  of  62  was  used  and  each 
measurement  elevation  was  corrected. 

Paltridge  and  Platt  (1976)  developed  the  Rayleigh  optical  depth  calculation 

tr{\,  h)  =  0.0088A("‘^'^^+°-2^)[exp(-0.1188/i  -  0.00116^^)]  (5.3) 

where  h  is  in  kilometers  and  A  is  in  microns.  Once  the  Rayleigh  optical  depth  is  known,  it 
is  assumed  that  haze  makes  up  the  difference  between  the  calculated  and  measured  optical 
depths. 

A  representative  surface  albedo  was  determined  from  Rowe  (1993)  and  the  solar 
spectred  flux  was  obtained  from  Iqbal  (1983). 

Figure  5.9  and  Figure  5.10  compare  the  measured  values  of  sky  light  radiance  and 
polarization  to  the  modeled  values  for  520  nm  and  670  nm.  The  radiance  values  from 
the  SSP  measurements  were  scaled  to  the  modeled  values  at  130°  elevation.  The  plots 
illustrate  the  model  in  three  stages  of  complexity:  first  with  only  Rayleigh  scattering  and 
no  surface  reflection,  the  second  step  adds  in  surface  albedo,  and  the  third  stage  has  both 
Rayleigh  and  haze  scattering  with  surface  reflection.  The  intensity  profile  requires  the 
addition  of  haze  into  the  model  before  a  realistic  prediction  is  achieved.  The  polarization 
is  reasonably  predicted  by  the  Rayleigh  scattering  case  with  surface  albedo.  The  addition 
of  haze  improves  the  model  prediction  at  all  angles.  Agreement  between  the  model  and 
measurements  is  generally  good,  except  in  the  vicinity  of  the  sun  where  values  of  negative 
polarization  were  observed.  Only  for  optically  thick  Rayleigh  atmospheres  did  the  model 
produce  negative  polarization  in  this  location.  Since  forward  scattering  of  the  modeled 
haze  is  unpolarized,  adding  haze  to  the  atmosphere  depolarized  the  negative  polarization 
induced  by  Rayleigh  scattering  near  the  sun.  Of  note  is  the  smaller  error  bar  for  130°  (90° 
above  the  solar  position)  when  compared  to  the  measurements  for  125  and  135°.  This  wzis 
a  direct  result  of  the  number  of  scans  taken  for  the  readings.  The  measurments  at  125 
and  135°  had  just  over  100  scans  of  data  and  the  130°  measurement  had  over  200  scans, 
thus  reducing  the  uncertainty.  Also,  the  maximum  errors  in  the  polarization  occur  close 
to  the  sun  and  the  horizon.  At  both  locations,  sky  light  intensity  changes  rapidly  with 
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viewing  angle.  The  most  probable  cause  for  this  error  is  SSP  channels  2  and  4  looking  in 
slighty  different  directions,  i.e.  the  parallel  channel  (2)  looking  above  the  perpendicular 
channel  (4). 


Number  of  Occurrences 
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Figure  5.8:  Radiance  values  error  distribution  for  130°  elevation.  Radiance  values  are 
grouped  into  20  equal  bins  between  the  maximum  and  minimum  readings  for  each  channel. 


Polarization  Radiance 


Polarization  Radiance 


Figure  5.10:  SSP  measured  intensity  and  polarization  compared  to  the  model  predicted 
values  for  670  nm. 


Chapter  6 


SUMMARY  AND  CONCLUSIONS 

The  SSP  is  in  its  third  iteration,  the  current  instrument  description,  optics  char¬ 
acteristics,  and  data  acquisition  systems  were  presented.  The  new  version  has  several 
improvements  over  the  previous  instruments: 

1.  Compact,  small  enough  to  fly  in  a  unmanned  aerospace  vehicle  (UAV). 

2.  Adds  polarization  capability.  This  version  has  the  capability  of  4  polarization 
channels;  parallel,  perpendicular,  and  right  and  left  hand  circular.  Only  the  parallel  and 
perpendicular  optics  were  installed  for  this  thesis  work.  Previous  versions  only  had  diffuse 
and  radiance  channels. 

3.  Increased  spectral  range.  In  the  current  configuration  the  SSP  spectral  range  is 
only  400  -  1100  nm.  When  the  new  detectors  are  installed  the  range  will  be  increased  to 
400  -  4000nm.  The  spectral  range  of  version  1  was  400  -  1200nm  and  version  2  was  400  - 
2500  nm 

The  SSP  is  ready  to  collect  atmospheric  data  from  airborne  platforms  at  this  time. 
The  instrument  hats  flown  test  missions  in  manned  aircraft  and  in  the  General  Atomics 
(GNAT)  UAV.  It  is  currently  in  compatibility  testing  for  the  Perseus-B  UAV  payload.  The 
SSP  will  provide  spectral  measurements  that  are  unique  to  the  remote  sensing  of  clouds, 
aerosols,  and  land  surfaces. 

Initial  calibration  wa^  conducted  here  at  Colorado  State  University.  The  calibration 
procedures,  to  include  the  optics  installation  and  alignment,  measuring  the  instruments 
field  of  view,  and  determining  the  angular  response  for  the  diffuse  channel  were  outlined. 
To  both  confirm  the  accuracy  of  the  calibration  and  a.s  a  study  of  the  sky  light  radiance, 
the  sky  light  intensity  and  polarization  were  measured  with  the  SSP  in  an  experiment  for 
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this  thesis.  A  polarized  radiative  transfer  model  was  utilized  to  examine  the  effects  local 
conditions  can  have  on  the  sky  light  intensity  and  polarization  and  to  evaluate  the  SSP 
data. 

6.1  Summary  of  Results 

6.1.1  Instrument 

The  optics  of  the  Scanning  Spectral  Polarimeter  were  characterized  in  this  paper. 

1.  Perpendicular  and  parallel  polarization  channels  are  aligned  90^  ±  .183®. 

2.  The  viewing  angle  for  the  three  narrow  field  of  view  channels  was  ^^20  mRad  for 
channels  2  and  3  and  ^19  mRad  for  channel  4.  These  values  are  slightly  lower  than  those 
predicted  by  ray  tracing. 

3.  The  angular  response  of  the  flashed  opal  diffuser  was  measured  as  a  function  of 
wavelength.  A  very  strong  forward  transmission  lobe  was  quantified  in  the  NIR  wave- 
lengths.  This  lobe  is  not  well  documented  by  the  manufacturer. 

4.  Calibration  coefficients  were  determined  and  allow  accurate  relationships  between 
spectral  radiances  measured  with  the  different  channels. 

6.1.2  Model 

The  polarized  radiative  transfer  model  was  used  to  isolate  effects  of  changing  Rayleigh 
optical  depth,  solar  position,  surface  reflection,  and  haze  loading.  Single  scatter  Rayleigh 
conditions  will  produce  an  atmosphere  with  0%  polarization  in  the  vicinity  of  the  sun  and 
100%  polarization  90®  from  the  sun. 

1.  As  Rayleigh  optical  depth  increases,  the  polarization  is  decreased  at  all  angles, 
and  can  even  become  negative.  Optically  thick  Rayleigh  atmospheres  have  a  significant 
area  of  negative  polarization  near  the  sun. 

2.  Increasing  surface  albedo  by  5%  causes  approximately  a  5%  decrease  in  the  maxi¬ 
mum  polarization  90®  from  the  sun.  The  effect  near  the  sun  is  much  smaller. 

3.  Increasing  the  solar  elevation  angle  causes  the  maximum  and  minimum  polariza¬ 
tions  to  shift  across  the  sky  and  also  a  slight  depolarization  at  all  angles. 
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4.  The  addition  of  haze  to  the  atmosphere  induces  strong  depolarization  at  all  levels. 
Haze  layers  of  relatively  small  optical  depth  (r/j  <  .1)  also  causes  the  intensity  of  diffuse 
radiation  to  increase  dramatically,  especially  near  the  sun. 

6.1.3  Experiment 

Polarization  measurements  made  with  the  SSP  were  in  generally  good  agreement  with 
values  predicted  by  the  model.  Rayleigh  scattering  alone  captures  the  major  elements  of 
the  sky  polarization  pattern.  Adding  surface  reflection  and  haze  layers  to  the  model  brings 
the  results  closer  to  that  measured  with  the  SSP.  The  average  polarization  error  in  the 
SSP  data  is  less  than  10%  with  respect  to  the  model.  Two  areas  of  maximum  disagreement 
were  close  to  the  sun,  where  the  model  predicted  0%  polarization  but  the  SSP  measured 
negative  polarization,  and  at  the  horizon.  The  maximum  error  for  both  cases  is  ~20%. 

6.2  Conclusions 

1.  The  SSP  provides  detailed  spectral  polarization  characteristics  suitable  for  research 
in  remote  sensing  of  the  earth’s  surface  and  the  atmosphere,  to  include  clouds  and  aerosols. 

2.  The  SSP  is  capable  of  measuring  the  spectral  sky  polarization  to  within  10% 
accuracy,  with  respect  to  model  predictions. 

3.  The  instrument  is  capable  of  resolving  the  effects  that  aerosol  scattering  and  surface 
reflection  have  on  sky  polarization.  The  non-spherical  properties  of  atmospheric  aerosols 
limit  the  degree  of  accuracy  to  which  the  model  can  predict  sky  polarization  patterns. 

4.  Characterization  of  the  flashed  opal  diffuser  spectral  transmission  function  will 
allow  the  SSP  to  determine  spectral  hemispheric  flux  within  5%  accuracy. 

6.3  Instrument  Improvements 

Several  instrument  improvements  are  currently  being  undertaken,  including: 

1.  The  NIR  detectors  are  being  replaced.  This  will  allow  measurements  from  1100  nm 
through  4000  nm. 

2.  The  optical  encoder  is  being  changed  so  measurements  are  recorded  at  even  in¬ 
tervals  around  the  CVF.  This  will  allow  all  detectors  to  make  measurement  at  identical 
wavelengths  and  simplify  polarization  calculations. 
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3.  The  detector  electronics  axe  being  upgraded  to  remove  the  DC  voltage  offset  prior 
to  signEil  amplification.  This  will  enable  the  data  part  of  the  signal  to  be  amplified  more, 
thus  reducing  the  significance  of  instrument  noise  on  errors. 


Bibliography 


Beiying,  W.  and  D.  Lu,  1988:  Retrieval  of  stratospheric  background  aerosol  scattering 
coefficient  from  twilight  polarization  measurements.  AppL  Optics,  27,  4899-4906. 

Deliver,  C.,  1988:  Luminance  and  polarization  of  the  sky  light  at  Seville  (Spain)  measured 
in  the  white  light.  Atmos,  Environ.,  22,  595-599. 

Chandrasekhar,  S.,  1960:  Radiative  Transfer.  Dover,  393  pp. 

Chandrasekhar,  S,  and  D.  D.  Ebert,  1954:  The  illumination  and  polarization  of  the  sunlit 
sky  on  Rayleigh  scattering.  Trans.  Am.  Philos.  Soc.,  44,  643-728. 

Coulson,  K.  L.,  1977:  Atmospheric  turbidity  determinations  by  skylight  measurements  at 
the  Mauna  Loa  Observatory.  Contrib.  Atmos.  Sci.,  13. 

Coulson,  K.  L.,  1980:  Characteristics  of  skylight  at  the  zenith  during  twilight  as  indicators 
of  atmospheric  turbidity.  1:  Degree  of  polarization.  Appl.  Optics,  19,  3469-3480. 

Coulson,  K.  L.,  1981:  Characteristics  of  skylight  at  the  zenith  during  twilight  as  indicators 
of  atmospheric  turbidity.  2:  Intensity  and  color  ratio.  Appl.  Optics,  20,  1516-1524. 

Coulson,  K.  L.,  1983:  Effects  of  the  El  Chichon  volcanic  cloud  in  the  stratosphere  on  the 
polarization  on  light  from  the  sky.  Appl.  Optics,  22,  1036-1050. 

Coulson,  K.  L.,  1988:  Polarization  and  Intensity  of  Light  in  the  Atmosphere.  A.  DEEPAK 
Publishing,  596  pp. 

Coulson,  K.  L.,  J.  V.  Dave,  and  Z.  Sekera,  1960:  Tables  Related  to  Radiation  Emerging 
from  a  Planetary  Atmosphere  with  Rayleigh  Scattering.  Univ.  of  California  Press, 
548  pp. 


88 


Dave,  J.  V.,  1970:  Intensity  and  polarization  of  the  radiation  emerging  from  a  plane- 
parallel  atmosphere  containing  monodispersed  aerosols.  Appl.  Optics^  9,  2673-2684. 

Dave,  J.  V.,  1978:  Review  paper:  Extensive  datasets  of  the  diffuse  radiation  in  realistic 
atmospheric  models  with  aerosols  and  common  absorbing  gases.  SoL  Energy^  21, 
361-369. 

Deirmendjian,  D.,  1969:  Electromagnetic  Scattering  on  Spherical  Poly  dispersions,  Amer¬ 
ican  Elsevier  Publishing  Co.,  290  pp. 

Deschamps,  P.,  F.  Breon,  M.  Leroy,  and  A.  Podaire,  1994:  The  POLDER  mission:  In¬ 
strument  characteristics  and  scientific  objectives.  IEEE  Trans.  Geosci.  Remote  Sens.,, 
32,  598-615. 

Deuze,  J.  L.,  F.  M.  Breon,  P.  Y.  Deschamps,  and  c.  Devaux,  1993:  Analysis  of  the 
POLDER  (POLarization  and  Directionality  of  Earth’s  Reflectances)  airborne  instru¬ 
ment  observations  over  land  surfaces.  Remote  Sens.  Environ..,  45,  137-154. 

Evans,  K.  F.  and  G.  L.  Stephens,  1991:  A  new  polarized  atmospheric  radiative  transfer 
model.  J.  Quant.  Spectosc.  Radiat.  Transfer,,  46,  413-423. 

Fitch,  B.  W.  and  K.  L.  Coulson,  1983:  Effects  of  skylight  at  south  pole  station,  Antarctica, 
by  ice  crystal  precipitation  in  the  atmosphere.  Appl.  Optics,,  22,  71-74. 

Goloub,  R,  J.  L.  Dueze,  M.  Herman,  and  Y.  Fouquart,  1994:  Analysis  of  the  POLDER 
polarization  measurements  performed  over  cloud  covers.  IEEE  Trans.  Geosci.  Remote 
Sens.,  32,  78-88. 

Grant,  1.  P.  and  G.  E.  Hunt,  1969:  Discrete  space  theory  of  radiative  transfer.  Proc.  Roy. 
Soc.  London,  A313,  183-197. 

Hitzfelder,  S.  J.,  G.  N.  Plass,  and  G.  W.  Kattawar,  1976:  Radiation  in  the  earth’s 
atmosphere:  Its  radiance,  polarization  and  ellipticity.  Appl.  Optics,  15,  2489-2500. 

Hovenier,  J.  W.,  1969:  Symmetry  relationships  for  scattering  of  polarized  light  in  a  slab 
of  randomly  oriented  particles.  J.  Atmos.  Sci.,  26,  488-499. 


89 


Iqbal,  M.,  1983:  An  Introduction  to  Solar  Radiation,  Academic  Press. 

Ishimaru,  A.,  D.  Lesselier,  and  C.  Yeh,  1984:  Multiple  scattering  calculations  for  non- 
spherical  particles  based  on  the  vector  radiative  transfer  theory.  Radio  Sci.^  19, 
1356-1366. 

Kattawar,  G.  W.,  G.  N.  Pllas,  and  S.  J.  Hitzfelder,  1976:  Multiple  scattered  radiation 
emerging  from  Rayleigh  and  continental  haze  layers.  1:  Radiance,  polarization,  and 
neutral  points.  Appl  Optics,,  15,  633-647. 

Meeus,  J.,  1991:  Astronomical  Algorithms,  Willmann-Bell  Inc,  429  pp. 

Nakajima,  T.  and  M.  D.  King,  1990:  Determination  of  the  optical  thickness  and  effective 
particle  radius  of  clouds  from  reflected  solar  radiation  measurements.  Part  I:  Theory. 
J.  Atmos.  Sci,.,  47,  1878-1893. 

Paltridge,  G.  W.  and  C.  M.  R.  Platt,  1976:  Radiative  Processes  in  Meteorology  and 
Climatology,  Elsevier  Sci.  Pub.  Co.,  318  pp. 

Rowe,  C.  M.,  1993:  Global  land-surface  albedo  modeling.  Int.  J.  Climat.^  13,  473-495. 

Scott,  J.  C.  and  G.  L.  Stephens,  1985:  A  visible-infrared  spectroradiometer  for  cloud 
reflectance  mea,surements.  J,  Phys.  E:  Sci.  Instrum..,  18,  697-701. 

Stackhouse,  P.  W.  and  G.  L.  Stephens,  1991:  A  theoretical  and  observational  study  of  the 
radiative  properties  of  cirrus:  Results  from  FIRE  1986.  J.  Atmos.  Sci,^  48,  2044-2059. 

Stackhouse,  P.  W.,  G.  L.  Stephens,  A.  J.  Heymsfield,  and  T.  Uttal,  1994:  Investigation  of 
the  effects  of  the  macrophysical  and  microphysical  properties  of  cirrus  clouds  on  the 
retrieval  of  optical  properties.  J.  Atmos.  Sci.  Submitted. 

Stephens,  G.  L.  and  J.  C.  Scott,  1985:  A  high  speed  spectrally  scanning  radiometer 
(SPERAD)  for  airborne  measurements  of  cloud  optical  properties.  J,  Atmos.  Ocean. 
Tech,,  2,  148-156. 


90 


Takano,  Y.  and  K.  Lion,  1989;  Solar  radiative  transfer  in  cirrus  clouds.  Part  I:  Single¬ 
scattering  and  optical  properties  of  hexagonal  ice  crystals.  J.  Atmos.  Sci.,  46,  3-22. 

Tsay,  S.,  P.  M.  Gabriel,  M.  D.  King,  and  G.  L.  Stephens,  1994:  A  fourier-riccati  approach 
to  radiative  transfer.  Part  II:  applications  to  spectral  reflectance  and  heating/cooling 
rates  in  the  5  December  1991  FIRE  IFO  cirrus  case  study.  J.  Atmos.  Sci.  Submitted. 

Wiscombe,  W.  J.,  1976:  On  initialization,  error  and  flux  conservation  in  doubling  method. 
J.  Quant.  Spectosc.  Radiat.  Transfer,  16,  477-489. 


Appendix  A 


CVF  CENTRAL  WAVELENGTHS 


Table  A.l:  Listing  of  the  central  wavelength  for  each  reading  made  by  the  SSP  as  a 
function  of  both  index  hole  number  and  degrees  of  rotation  of  the  CVF  wheel. 


Detector  Angle 

240 

120 

60 

180 

0 

Hole 

Angle 

Det  1 

Det  2 

Det  3 

Det  4 
_ 

Det  5 

Det  6 

1 

3 

2727.4 

1831.67 

842 

509 

1355 

4026.07 

2 

6 

2662.47 

1795 

814.5 

521.5 

1330 

3961.13 

3 

9 

2597.53 

1758.33 

788.5 

535 

1305 

3896.2 

4 

12 

2532.6 

1721.67 

761.5 

548 

1280 

3831.27 

5 

15 

2467.67 

1685 

737.5 

561 

1255 

3766.33 

6 

18 

2402.73 

1648.33 

712.5 

575 

1230 

3701.4 

7 

21 

2337.8 

1611.67 

687 

588 

1204 

3636.47 

8 

24 

2272.87 

1575 

664 

601 

1179 

3571.53 

9 

27 

2207.93 

1538.33 

639 

614 

1154 

3506.6 

10 

30 

clear 

1501.67 

clear 

627 

1129 

3441.67 

11 

33 

2198.33 

1465 

381 

640 

1104 

3376.73 

12 

36 

2161.67 

1428.33 

394 

653.5 

1077 

3311.8 

13 

39 

2125 

1391.67 

407 

666.5 

1050 

3246.87 

14 

42 

2088.33 

1355 

420 

680 

1025 

3181.93 

15 

45 

2051.67 

1318.33 

433 

693 

996 

3117 

16 

48 

2015 

1281.67 

445 

706 

970 

3052.07 

17 

51 

1978.33 

1245 

458 

719.5 

945 

2987.13 

18 

54 

1941.67 

1208.33 

470.5 

733 

919 

2922.2 

19 

57 

1905 

1171.67 

483. 

742 

892 

2857.27 

20 

60 

1868.33 

clear 

496 

clear 

867 

2792.33 

21 

63 

1831.67 

1355 

509 

4026.07 

842 

2727.4 

22 

66 

1795 

1330 

521.5 

3961.13 

814.5 

2662.47 

23 

69 

1758.33 

1305 

535 

3896.2 

788.5 

2597.53 

24 

72 

1721.67 

1280 

548 

3831.27 

761.5 

2532.6 

25 

75 

1685 

1255 

561 

3766.33 

737.5 

2467.67 

26 

78 

1648.33 

1230 

575 

3701.4 

712.5 

2402.73 

27 

81 

1611.67 

1204 

588 

3636.47 

687 

2337.8 

28 

84 

1575 

1179 

601 

3571.53 

664 

2272.87 

92 


93 


Table  A.l  (cont) 


Detector  Angle 

300 

240 

120 

60 

180 

0 

Hole 

Angle 

Det  1 

Det  2 

Det  3 

Det  4 

Det  5 

Det  6 

69 

207 

639 

614 

2207.93 

1538.33 

3506.6 

1154 

70 

210 

clear 

627 

clear 

1501.67 

3441.67 

1129 

71 

213 

381 

640 

2198.33 

1465 

3376.73 

1104 

72 

216 

394 

653.5 

2161.67 

1428.33 

3311.8 

1077 

73 

219 

407 

666.5 

2125 

1391.67 

3246.87 

1050 

74 
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420 
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1355 

3181.93 

1025 

75 
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433 

693 

2051.67 

1318.33 

3117 

996 

76 
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1281.67 

3052.07 
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77 
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2987.13 

945 

78 
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470.5 

733 

1941.67 

1208.33 

2922.2 

919 

79 

237 

483 

742 

1905 

1171.67 

2857.27 

892 

80 

240 

496 

clear 

1868.33 

clear 

2792.33 

867 

81 

243 

509 

4026.07 

1831.67 

1355 

2727.4 

842 

82 
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521.5 

3961.13 

1795 

1330 

2662.47 

814.5 

83 

249 

535 

3896.2 

1758.33 

1305 

2597.53 

788.5 

84 
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548 

3831.27 

1721.67 

1280 

2532.6 

761.5 

85 

255 

561 

3766.33 

1685 

1255 

2467.67 

737.5 

86 

258 

575 

3701.4 

1648.33 

1230 

2402.73 

712.5 

87 

261 

588 

3636.47 

1611.67 

1205 

2337.8 

687 

88 
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3571.53 

1575 

1180 

2272.87 
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3506.6 

1538.33 

1154 

2207.93 
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90 
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3441.67 

1501.67 

1128 
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91 

275.3 

650 

3326.95 

1436.89 

1082 

2170.22 

395.38 

92 

280.6 

673 

3212.24 

1372.11 

1037 

2105.44 

417.76 

93 

285.9 

696 

3097.52 

1307.33 

986 

2040.67 

440.13 

94 

291.2 

718 

2982.8 

1242.56 

943 

1975.89 

462.51 

95 

296.5 

740 

2868.09 

1177.78 

896 

1911.11 

484.89 

96 

301.8 
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2753.37 

clear 

848.5 

1846.33 

507.27 

97 

307.1 

3937.32 

2638.66 

1320 

803.5 

1781.56 

529.64 

98 

312.4 

3822.61 

2523.94 

1275 

758 

1716.78 

552.02 

99 

317.7 

3707.89 

2409.23 

1230 

714 

1652 

574.40 

100 

323 

3593.18 

2294.51 

1182 
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1587.22 

596.78 

101 

328.3 

3478.46 
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1143 

clear 

1522.44 

619.16 

102 

333.6 

3363.75 

2191 

1098 

390 

1457.67 

641.53 

103 

338.9 

3249.03 

2126.22 

1052 

411 

1392.89 

663.91 

104 

344.2 

3134.32 

2061.44 

1008 

432.5 

1328.11 

686.29 

105 

349.5 

3019.6 

1996.67 

958 

453.5 

1263.33 

708.67 

106 

354.8 

2904.88 

1931.89 

912 

475 

1198.56 

731.04 

107 

360 

2792.33 

1868.33 

867 

496 
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PROCEDURES  FOR  CALIBRATION  EXERCISES 
B*1  Optics  Alignment 

The  Glan-Taylor  polarization  cube  in  the  perpendicular  channel  must  be  aligned  90° 
from  the  Glan-Taylor  in  the  parallel  channel.  Equipment  layout  for  aligning  the  optics  is 
illustrated  in  Figure  3.1.  Optical  channels  2  and  4  are  used  for  the  polarization  channels 
since  they  are  physically  180°  apart  from  each  other  in  the  optics  tower. 

B.1.1  Alignment  Procedure 

1.  Remove  the  optical  tower  from  the  SSP  vacuum  chamber  and  mount  on  the  jig 
with  channel  2  located  in  the  top  right  position. 

2.  Remove  both  the  focusing  lens  (from  optics  layer  3)  and  Glan-Taylor  polarizing 
cubes  (from  optics  layer  1)  if  they  are  installed  for  optical  channels  2  and  4. 

3.  Put  a  centering  guide  in  the  back  (layer  3)  of  the  optics  tower.  Align  the  laser 
beam  on  the  center  of  the  guide. 

4.  Put  a  centering  guide  in  the  front  (layer  1)  of  the  optics  tower.  The  laser  beam 
should  hit  this  new  guide  in  its  center.  If  the  laser  beam  does  not  hit  the  second  guide 
in  its  center,  the  optics  tower/laser  beam  are  not  aligned  with  each  other.  Adjust  the 
position  of  the  optics  tower  for  large  corrections.  Make  small  corrections  with  the  two 
mirrors.  Use  mirror  1  to  adjust  the  originating  position  of  the  laser  beam  on  mirror  2 
and  mirror  2  to  adjust  the  beam  direction.  Repeat  steps  3  and  4  until  the  laser  beam  is 
centered  on  both  the  back  and  front  guides. 

5.  After  the  beam  is  aligned  on  both  centering  guides,  install  the  iris.  Close  the 
aperture  on  the  iris  so  the  opening  is  slightly  larger  than  the  laser  beam.  Verify  alignment 
of  the  optical  tower  by  installing  the  Glan-Taylor  polarizing  cube.  The  Glan-Taylor  will 
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reflect  some  of  the  incident  light.  Look  on  the  back  of  the  iris  to  see  the  direction  of 
the  reflected  beam.  The  reflected  beam  should  be  cJigned  with  the  incident  beam.  Small 
imperfections  in  the  optics  tower  or  lens  construction  will  cause  this  alignment  to  be 
slightly  off.  This  error  w2ls  measured  to  be  0.28°  between  the  incident  and  reflected  beam 
for  channel  2,  and  approximately  half  this  value  for  channel  4.  Once  the  laser  beam  is 
aligned  with  the  optical  tower,  do  not  adjust  the  laser  beam  or  the  optical  tower  jig.  The 
following  adjustments  will  be  accomplished  by  moving  the  detector  or  lens  holders. 

6.  Remove  the  Glan-Taylor  and  the  centering  guides  from  the  optics  tower.  Set  up 
the  detector,  centering  it  on  the  laser  beam.  Also  ensure  the  light  reflected  by  the  detector 
is  aligned  with  the  incident  beam.  Maximum  signal  from  the  detector  well  be  obtained 
when  the  detector  is  properly  centered  and  aligned. 

7.  Put  a  lens  centering  guide  in  lens  holder  2.  Center  the  guide  in  the  laser  beam. 
Once  the  lens  is  centered  it  can  be  aligned  by  putting  a  piece  of  paper  behind  the  guide.  If 
aligned  properly,  a  circle  of  light  will  be  transmitted.  If  the  guide  is  out  of  alignment  it  will 
create  an  interference  pattern.  To  align  the  guide,  remove  the  interference.  This  will  yield 
only  a  crude  alignment.  Remove  the  centering  guide  and  install  the  Glan-Taylor.  Fine 
tune  the  alignment,  by  adjusting  the  reflected  beam  so  it  is  coincident  with  the  incident 
beam.  Remove  the  Glan-Taylor  and  install  the  centering  guide  to  verify  the  lens  holder  is 
still  centered  on  the  laser  beam.  Repeat  the  process  until  the  lens  holder  is  both  centered 
and  aligned. 

8.  Repeat  step  7  for  lens  holder  1. 

Once  all  the  equipment  is  aligned  with  the  laser  beam,  Glan-Taylor  polarization  cubes 
are  ready  for  installation  into  the  SSP  optical  tower.  There  are  several  methods  to  align 
the  Gian- Taylors.  These  methods  normally  include  setting  a  Glan-Taylor  at  a  reference 
angle,  shooting  the  laser  beam  through  this  reference  cube  and  then  through  the  Glan- 
Taylor  being  aligned.  Either  the  maximum  signal  (the  two  cubes  at  the  same  angle)  or  the 
minimum  signal  (the  two  cubes  90°  apart)  are  sought.  For  the  equipment  setup  used  it 
was  very  difficult  to  determine  the  maximum  signal  since  the  detector  became  saturated 
and  readings  near  the  maximum  were  very  unstable.  In  an  attempt  to  overcome  this 
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problem,  a  filter  was  used  to  attenuate  the  signal,  but  readings  from  the  detector  were  not 
reliable.  To  get  a  discernible  change  in  the  signal,  the  polarization  cube  had  to  be  rotated 
3-5®.  On  the  other  hand,  detecting  the  minimum  signal  was  very  accurate.  The  base-line 
signal  from  the  detector  with  no  light  hitting  the  detector  was  0.48mV.  The  detector  was 
sensitive  to  any  incoming  light.  During  this  phase  of  the  calibration  all  external  light  was 
eliminated.  The  lights  were  off,  windows  were  covered,  and  the  doors  were  shut  with  the 
hole  at  the  bottom  covered.  The  sensitivity  at  this  range  was  measured  at  more  than 
lOOmV  per  degree  of  rotation.  The  following  alignment  procedure  utilizes  this  extreme 
sensitivity  at  determining  the  minimum  signal  when  the  two  Glan-Taylors  are  at  90°  to 
each  other. 

9.  Insert  a  Glan-Taylor  into  lens  holder  1.  This  is  the  reference  cube.  This  cube  must 
be  set  square  with  the  light  table.  It  was  aligned  to  allow  radiation  perpendicular  to  the 
light  table  to  pass. 

10.  Insert  a  Glan-Taylor  into  channel  2  of  the  optics  tower.  Rotate  the  Glan-Taylor 
until  the  minimum  signal  is  obtained.  Lock  the  Glan-Taylor  into  position.  This  is  difficult, 
as  even  a  slight  rotation  of  the  Glan-Taylor  caused  by  any  vibration  or  uneven  tightening 
of  the  lens  clamp  will  change  the  lens  orientation  and  thus  cause  the  transmitted  signal 
to  vary.  Once  the  minimum  signal  is  obtained,  the  Glan-Taylor  in  the  optics  tower  will 
be  aligned  90°  from  the  reference  cube.  Channel  2  will  now  eliminate  all  perpendicular 
radiation  and  therefore  detect  parallel  polarization.  The  minimum  reading  achieved  was 
3.16mV. 

11.  Rotate  the  optics  tower  180°.  This  will  move  channel  4  into  the  laser  beam  path. 
Verify  alignment  of  the  optics  tower  and  laser  beam  by  placing  the  centering  guides  and 
Glan-Taylors  into  optical  channel  4.  No  adjustments  should  be  required.  If  the  optics 
tower/laser  beam  are  not  aligned  after  the  rotation,  the  optics  jig  has  most  likely  moved 
and  the  entire  process  must  be  repeated. 

12.  Remove  the  centering  guides  and  Glan-Taylor  from  optical  channel  4.  Insert  a 
Glan-Taylor  into  lens  holder  2.  Rotate  this  Glan-Taylor  until  the  minimum  transmitted 
signal  is  obtained.  It  is  easier  to  adjust  this  lens  since  the  lens  holder  makes  it  simpler 
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to  lock  the  lens  into  position.  Once  the  minimum  signal  is  obtained,  the  Glan-Taylor  in 
lens  holder  2  will  be  aligned  90°  from  the  reference  Glan-Taylor  in  lens  holder  1  and  thus 
will  allow  radiation  parallel  to  the  light  table  to  pass.  The  minimum  signal  obtained  was 
O.TOmV.  The  Glan-Taylor  in  lens  holder  2  now  becomes  the  reference  polarization  cube. 

13.  Remove  the  Glan-Taylor  from  lens  holder  1  and  install  it  in  the  optics  tower,  layer 
1  of  channel  4.  Repeat  the  alignment  procedure  of  step  10.  The  Glan-Taylor  in  channel  4 
will  now  be  aligned  90°  from  the  one  in  channel  2.  Channel  4  will  detect  perpendicular 
polarization.  The  minimum  signal  obtained  was  2.07mV. 

B.1.2  Errors  in  alignment 

To  calculate  the  error  in  aligning  the  two  Glan-Taylors,  the  voltages  from  the  three 
settings  were  added  together  and  three  times  the  base-line  voltage  was  then  subtracted. 
This  represents  the  maximum  error  since  it  assumes  all  errors  in  setting  the  lenses  were 
in  the  same  direction.  This  error  is  4.49mV.  Using  the  sensitivity  of  lOOmV  per  degree, 
measured  earlier,  the  error  is  .045°.  Also,  while  re-positioning  the  optical  tower,  it  was 
noted  the  jig  allowed  the  optical  tower  to  rotate  slightly.  This  movement  was  measured 
to  be  0.138°.  This  error  is  then  added  to  the  error  calculated  above  providing  the  final 
positioning  of  the  two  Glan-Taylors  as  90°  ±  .183°.  This  is  a  conservative  estimate,  since 
all  errors  are  assumed  to  be  in  the  same  direction  and  as  the  optics  tower  is  tightened  to 
the  jig,  movement  of  the  optics  tower  is  reduced. 

B.2  Field  of  View 

The  purpose  of  this  exercise  is  to  determine  the  field  of  view  (FOV)  for  the  three 
narrow  field  of  view  (NFOV)  channels  (channels  2-4).  A  plan  view  of  the  equipment 
setup  used  for  this  exercise  is  shown  in  Figure  3.3  and  the  side  view  is  shown  in  Figure 
B.l.  The  SSP  is  set  on  a  leveling  bench  and  the  light  source  and  laser  sit  on  scissors  jacks. 
The  light  source  is  on  a  rolling  bench.  This  exercise  is  conducted  in  a  long  hallway. 

B,2.1  FOV  Measurement  Procedure 


1.  Center  the  table  in  the  hallway.  Level  the  table  and  then  the  leveling  bench. 
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Figure  B.l:  FOV  measurement  layout. 

2.  Set  laser  and  SSP  on  the  leveling  bench,  SSP  in  the  back  and  Icuser  in  the  front. 
Align  the  SSP  center  line  in  the  vertical.  Adjust  the  height  of  the  laser  to  the  height  of 
the  entrance  window  for  channel  2. 

3.  Shoot  the  laser  down  the  center  of  the  hallway  towards  the  light  source.  Position 
the  light  source  in  the  center  of  the  hallway  and  adjust  the  height  of  the  light  source  so 
the  laser  beam  hits  it  in  the  center  of  the  element, 

4.  Move  the  SSP  into  position  on  the  leveling  bench.  Move  the  laser  to  the  light 
source  end  of  the  hallway,  behind  the  light  source. 

5.  Position  the  laser  in  its  new  position  so  it  is  in  on  the  experiment  center  line, 
shooting  down  the  center  of  the  hallway  towards  the  SSP.  Adjust  the  height  of  the  laser 
to  the  height  of  the  center  of  the  light  source  element.  Temporarily  lower  the  light  source. 

6.  Align  the  latser  so  the  laser  beam  hits  the  SSP  channel  2  entrance  window. 

7.  Align  the  SSP  on  the  leveling  table  so  the  laser  beam  is  reflected  back  along  its 
incident  path.  Lock  the  SSP  into  position  on  the  leveling  bench. 

8.  Adjust  the  light  source,  height  and  position,  so  the  laser  hits  the  center  of  the 
element.  The  line  from  the  SSP  to  the  light  source  now  defines  the  experiment  center  line. 
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9.  Mark  the  position  of  the  light  source  bench.  Tape  a  yard  stick  or  tape  measure  to 
the  floor  perpendicular  to  the  experiment  centerline.  This  will  be  used  to  determine  the 
distance  the  light  source  is  from  the  centerline. 

10.  Measure  the  distance  between  the  SSP  and  light  source,  this  distance  was  20  feet 
for  our  setup.  At  20  feet  we  received  ample  signal  to  conduct  the  experiment,  yet  not 
saturate  the  detector  and  there  was  enough  room  in  the  hallway  to  move  the  light  source 
out  of  the  field  of  view  of  the  SSP. 

11.  Attach  a  cardboard  screen  to  the  front  of  the  light  source  bench.  This  screen 
moves  with  the  light  source  and  will  stop  most  the  reflections  from  the  walls.  Cut  a  small 
hole  for  the  light  to  pass  through  and  reach  the  SSP.  Use  the  laser  to  find  the  center  of 
the  hole.  The  hole  must  be  large  enough  for  light  from  the  source  to  illuminate  the  SSP 
when  the  source  is  moved  to  the  side  of  the  hallway. 

12.  To  reduce  reflections  further,  hang  black  felt  curtains  across  the  hallway.  Cut  a 
hole  in  the  curtains  to  allow  light  from  the  source  to  reach  the  SSP. 

13.  Hang  a  black  felt  curtain  behind  the  laser  to  reduce  the  reflections  from  the  back 

wall. 

14.  Turn  off  the  laser. 

15.  The  light  source  is  now  moved  along  the  perpendicular  to  the  experiment  center 
line.  Light  intensity  received  at  the  SSP  is  measured  as  the  distance  from  the  center  line 
in  varied. 

16.  Repeat  for  each  narrow  field  of  view  channel. 

B.3  Cosine  Response  of  Diffuser 

The  object  of  this  exercise  is  to  measure  how  the  transmission  of  the  diffuser  lens  on 
channel  1  varies  as  the  angle  of  incident  light  changes.  Setup  is  similar  to  that  for  the 
narrow  field  of  view  mea.surement,  however,  the  distance  between  the  SSP  and  the  light 
source  is  only  53.3  cm.  A  plan  view  of  the  experiment  setup  is  shown  in  Figure  3.10.  The 
side  view  is  very  similar  to  that  shown  for  the  FOV  exercise  in  Figure  B.l,  except  the 
distances  are  shorter  and  the  curtains  are  positioned  differently  to  reduce  reflections. 
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B.3.1  Diffuser  Transmission  Function  Measurement  Procedure 

L  Level  table  and  leveling  bench.  Align  the  SSP  centerline  in  the  vertical. 

2.  Set  the  laser  on  the  leveling  bench  at  the  height  of  the  diffuser  (channel  1)  entrance 
lens.  Shoot  the  laser  towards  the  source  light.  Adjust  the  height  of  the  light  so  the  laser 
beam  hits  the  center  of  the  filament. 

3.  Set  the  SSP  on  the  leveling  bench  with  the  instrument  center  line  in  the  vertical. 
Draw  a  vertical  line  from  the  center  of  the  diffuser  lens  to  the  base  of  the  SSP  stand. 
This  line  will  help  keep  the  SSP  in  the  proper  position  as  it  is  rotated.  Align  the  SSP 
square  with  the  leveling  bench.  Draw  a  line  across  the  leveling  bench,  parallel  with  the 
SSP  stand.  Set  the  front  of  the  SSP  stand  on  this  new  line.  Mark  the  point  on  the  leveling 
bench  directly  below  the  center  of  the  diffuser  lens.  The  angle  of  SSP  rotation  will  be 
measured  with  respect  to  the  leveling  bench,  so  this  must  be  as  exact  as  possible. 

4.  Move  the  laser  behind  the  light  source.  Adjust  and  level  the  laser  beam  to  the 
height  of  the  center  of  the  light  filament.  Temporarily  lower  the  source  light. 

5.  Shoot  the  laser  towards  the  SSP  hitting  the  diffuser  lens  for  channel  1.  Adjust  the 
laser  so  the  reflected  beam  lines  up  with  the  incident  beam.  Do  not  move  the  SSP  for  this 
adjustment. 

6.  Move  the  light  source  into  the  path  of  the  laser  and  adjust  the  height  so  the  laser 
is  centered  on  the  filament. 

7.  Set  black  felt  curtains  to  reduce  reflection  from  the  walls. 

8.  The  SSP  head  is  rotated  from  0°  (looking  straight  at  the  light)  through  90°  (looking 
perpendicular  to  the  light).  Intensity  measurements  are  recorded  as  the  angle  is  varied. 
The  angle  of  the  SSP  was  measured  with  a  protractor.  Care  must  be  taken  so  the  distance 
from  the  light  source  and  the  position  of  the  center  of  the  diffuser  lens  do  not  change  as 
the  SSP  is  rotated.  Use  the  lines  drawn  in  step  3  to  keep  the  diffuser  lens  of  the  SSP  in 
the  proper  position  as  it  rotated. 

9.  Measurements  were  recorded  every  degree  from  0°  through  5°,  then  every  5° 
from  5°  through  80°.  Measurements  between  0°  are  5°  important  to  define  the  forward 
transmission  lobe  for  the  NIR  wavelengths. 
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B.4  Detector  Response,  Narrow  Field  of  View  Channels 

This  is  the  most  critical  part  of  the  calibration.  Equipment  setup  for  this  exercise 
is  illustrated  in  Figure  3.17.  This  exercise  was  conducted  in  a  large  room.  The  radiance 
measured  by  a  channel  must  be  related  to  the  radiance  measured  by  the  same  channel 
at  a  different  gain  setting  and  to  the  radiance  measured  by  the  other  channels  at  the 
same  and  different  gain  settings.  To  increase  the  accuracy  of  the  intensity  relationships  all 
three  NFOV  channels  detectors  were  calibrated  at  the  same  time.  This  introduces  slightly 
different  geometry  for  each  channel.  The  measurements  in  Figure  3.17  axe  shown  for  the 
center  of  the  SSP.  Exact  angles  and  distances  must  be  calculated  for  each  channel  of  the 
SSP. 

B.4.1  NFOV  Detector  Response  Measurement  Procedure 

1.  Level  table  and  leveling  bench. 

2.  Set  the  SSP  on  the  leveling  bench,  position  and  alignment  are  not  critical. 

3.  Set  the  laser  on  a  second  table.  Level  the  Iciser  and  adjust  its  height  to  the  height 
of  the  center  of  the  SSP. 

4.  Shoot  the  laser  towards  the  lambertian  plate.  The  lambertian  plate  worked  just  as 
advertised,  the  laser  bean  was  reflected  evenly  in  all  directions.  To  align  the  lambertian 
plate,  attach  a  plexiglass  triangle  to  its  front.  Put  a  piece  of  tape  on  the  back  of  the  triangle 
where  the  laser  bean  hits.  The  laser  beam  will  be  reflected  for  aligning  the  equipment. 

5.  Adjust  the  height  and  position  of  the  lambertian  plate  so  the  laser  beam  hits  the 
center  of  the  plate  and  then  reflects  from  the  plexiglass  triangle  towards  the  SSP. 

6.  Adjust  the  SSP  position  so  the  laser  beam  reflected  from  the  lambertian  plate  hits 
it  in  the  center.  Align  the  SSP,  by  putting  a  second  plexiglass  triangle  on  the  front  of  the 
SSP  and  reflect  the  laser  beam  along  its  incident  path. 

7.  Measure  the  angles  and  distances  between  the  center  of  the  lambertian  plate  and 
SSP.  Also  measure  the  angle  between  the  lambertian  plate  and  laser. 

8.  Adjust  the  source  lamp  height  and  position  so  the  laser  hits  the  center  of  the 


element. 
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9.  Measure  the  distance  between  the  light  source  and  the  center  of  the  lambertian 
plate. 

10.  Turn  off  the  laser. 

11.  Measure  the  signal  received  by  the  SSP  at  each  gain  setting  used  to  collect  solar 


data. 
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DIFFUSER  TRANSMISSION  CURVES 


Diffuser  transmission  curve  parameters  for  the  equation: 


Receiver  Power  (6) 
ReceivedPower{6  =  0) 


=  +  60  +  c  +  d{exp[—e6]) 


(C.l) 


Table  C.l:  Diffuser  transmission  coefficients  for  equation  C-1. 


Hole 

Wavelength  (nm) 

a 

b 

c 

d 

e 

81 

509.00 

-1.55E-04 

1.07E-03 

8.94E-01 

0 

0 

82 

521.50 

-1.43E-04 

-1.50E-03 

9.48E-01 

0 

0 

83 

535.00 

-1.25E-04 

-1.42E-03 

9.28E-01 

0 

0 

84 

548.00 

-1.40E-04 

-1.17E-03 

9.82E-01 

0 

0 

85 

561.00 

-1.33E-04 

-1.45E-03 

9.49E-01 

0 

0 

86 

575.00 

-1.35E-04 

-1.54E-03 

9.76E-01 

0 

0 

87 

588.00 

-1.34E-04 

-1.71E-03 

9.90E-01 

0 

0 

88 

601.00 

-1.33E-04 

-1.79E-03 

9.90E-01 

0 

0 

89 

614.00 

-1.32E-04 

-1.94E-03 

9.91E-01 

0 

0 

90 

627.00 

-1.37E-04 

-1.86E-03 

l.OOE+00 

0 

0 

91 

650.00 

-1.31E-04 

-2.21E-03 

9.99E-01 

0 

0 

92 

673.00 

-1.26E-04 

-2.79E-03 

l.OOE+00 

0 

0 

93 

696.00 

-1.20E-04 

-3.26E-03 

9.99E-01 

0 

0 

94 

718.00 

-1.18E-04 

-3.45E-03 

l.OOE+00 

0 

0 

95 

740.00 

-1.18E-04 

-3.53E-03 

9.86E-01 

0 

0 

68 

664.00 

-1.20E-04 

-2.94E-03 

9.97E-01 

0 

0 

67 

687.00 

-1.28E-04 

-2.94E-03 

l.OOE+00 

0 

0 

66 

712.50 

-1.23E-04 

-3.19E-03 

l.OlE+00 

0 

0 
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